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INTRODUCTION 


This  report  covers  the  work  accomplished  under  AFOSR 
grant  number  AFOSR  414-67  which  has  not  yet  been  covered  in 
our  existing  reports  or  publications,  in  particular,  we  have 
attempted  to  present  the  theory  and  empirical  results  which 
we  believe  represent  the  state-of-the-art  in  travel-time 
estimation  and  the  determination  of  epicenters.  Many  of  the 
findings  follow  from  results  previously  reported.  These 
results  are  not  discussed  in  detail  in  this  report;  thus,  the 
reader  may  need  to  refer  to  some  of  the  references  in  order  to 
follow  the  discussion.  Many  of  the  pertinent  reference  papers 
are  in  the  "Special  Number-1968  Seismological  Tables  for  P  Phases" 
Bulletin  of  bhe  Seismological  Society  of  America .  Vol.  58,  No.  4. 

In  the  report  we  first  develop  a  useful,  mathematical 
model  for  the  teleseismic  travel-times  and  then  discuss  the  sta¬ 
tistical  techniques  involved  in  the  estimation  of  model 
parameters.  The  problems  which  limit  the  accuracy  of  epicenter 
determinations  are  discussed  in  detail,  we  bslieve  that  the 
theory  and  numerical  results  presented  in  this  report  provide 
a  good  appraisal  of  current  capabilities  for  epicenter  determin¬ 
ation  and  a  satisfactory  base  for  further  studies  designed  to 
improve  the  accuracy  of  estimated  location  parameters. 


SEISMIC  TRAVEL-TIME  MODELS  FOR  TELES EISMIC  P-WAVES 

Basic  to  any  statistical  technique  is  the  model  used 
to  represent  the  travel  times.  This  model  must  adequately 
explain  differences  between  observed  and  predicted  arrival 
times,  in  order  to  be  useful  the  model  must  (1)  have  a 
reasonable  physical  basis  and  (2)  be  estimable  in  the  sta¬ 
tistical  sense.  The  first  requirement  enables  the  model  to 
be  related  to  underlying  theory  and  the  second  enables  per¬ 
tinent  parameters  to  be  estimated  from  observed  data,  we 
begin  with  a  complete  definition  of  such  a  model  and  then 
consider  approximations  needed  to  make  the  location  prob¬ 
lem  and  the  estimation  of  travel  times  tractable. 

Our  general  model  is: 

°j  +TC^j'hj)+^j  +  eAJk  (1) 

where  i  refers  to  the  station, 
j  refers  to  the  source, 
k  refers  to  the  observation, 
fcijk  the  observed  arrival  time, 

Oj  is  the  origin  time  of  the  source, 

T(^ij '  hj)  is  the  predicted  travel  time  from  the 
lib  source  to  the  ith  station  (from  double-entry 
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tables ) , 

A  *  U  the  source  to  station 

Un  great  Ci^cle  degrees), 

hj  ^  thS  S°UrCe  dePth  kilometers) , 

{*S)  ij  ^  M  ^nteracfcion  time  term,  and 

€ 

ijk  ls  the  error  term. 

In  this  model  € . . 

ex  .  1]k  ei'r0r  °f  Ovation  and  can  be 

examined  by  holding  i  and  j  fixed  and  1  v 

d  lo°kmg  at  repeated 
observations  (varying  k) .  Work  „ 

of  Freedman  (1968,  1966) 
indicates  that  i-hi  -  ' 

*  error  term  can  be  approximated  by  a 

mixture  of  normal  distributions.  The  •  . 

principal  error  com- 

nen  s  are.-  (Freedman,  1968) 

(1)  Miscounts  -  These  «« 

These  occur  m  multiples  of  seconds, 

minutes,  quarter  hours 

hours,  or  even  hours,  independent 

or  the  errors  listed  below. 

<2)  Mis identifications  -  These  arise  as 

iSe  as  a  result  of 

poor  signal-to-noise  ratios  mb,- 

•  This  error  and  the 

fading  error  are  dependent  variables. 

<3’  Instrumental  errors  -  These  arise  . 

arise  from  variations 

- 

ment  response,  etc. 
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(4)  Reading  errors  -  This  is  the  residual  error  which 
would  remain  even  though  the  above  were  eliminated. 

By  employing  only  the  better  stations  and  larger  events 
it  is  possible  to  virtually  eliminate  types  (2)  and  (3).  m 
practice,  it  can  then  be  assumed  that  any  £.jk  consists  of 
components  (1)  and  (4).  This  implies,  unfortunately,  that 
the  eijk  are  contaminated  (Tukey, 1960,  1962).  since  contam¬ 
ination,  even  a  very  small  proportion,  may  vitiate  the  usual 
estimators,  data  sets  must  routinely  be  truncated,  with  the 
data  truncated  for  outliers,  then  the  usual  least  squares 
techniques  are  applicable  (Tukey,  I960;  Dixon,  1953).  We 
will  return  to  this  point  later.  However,  with  truncation 
it  is  possible  to  assume  that: 

eijk“(Truncated)HID  (0,62)  over  all  i,  j,  k.  (2) 
Studies  employing  the  above  considerations  and  using  repeat¬ 
ed  explosions  at  the  Nevada  test  site,  indicate  that  ® 2  is 
of  the  order  of  0.01  to  0.04  sec2. 

Furthermore,  we  shall  place  restrictions  on  the  travel 
time  model  based  on  geophysical  considerations.  We  require, 

from  the  concept  of  reciprocity,  that  all  travel  time  terms 

be  symmetrical;  that  is, 

= 


(3) 


-4- 


The  travel  time  remains  the  same  if  we  interchange 
and  station.  In  addition,  we  restrict  the  value  o 
such  that, 


source 


Z  5  9  *  A; j  >  100°. 


The  model  given  by  equation  (1)  (Model  1)  meets  the 
two  criteria.  Consider  a  fixed  source,  j,  and  station,  i. 
The  observed  arrival  time  with  repeated  events  at  j  consists 
of  an  origin  time  plus  a  calculated  travel  time  (from  double 
entry  tables)  plus  a  correction  term  for  the  difference  in 
the  true  travel  time  from  j  to  i  and  the  tabled  value  and 
finally,  plus  measurement  error.  Also  with  explosion  data 
it  is  possible  to  estimate  the  only  important  unknown, 

(0  3)^.  With  the  assumption  of  normality  the  unique  mini¬ 
mum  variance  unbiased  estimator  of  (0(7)  .  .  is 

jtiUjk-Qjk  ~T(A;jjhj))/n 


where  i,  j  are  fixed, 

n  =  the  number  of  events  at  source  j,  and 
Ojk  is  the  origin  time  for  the  kth  event  at  j. 

Since  in  practice  it  is  difficult  to  repeat  explosions  at 
exactly  the  same  site  another  way  to  accomplish  the  estima¬ 
tion  is  to  have  one  event  and  a  number  of  recorders  at  each 
station.  The  minimum  variance  unbiased  estimator  is  again 
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given  by  equation  (4)  with  0jk  replaced  by  0j.  The  explo¬ 
sion  data  from  the  Nevada  test  site  could  reasonably  fit 
either  of  these  situations.  In  fact  the  procedure  amounts 
to  a  calibration  of  the  scurce-station  pair. 

In  practice,  the  procedure  cannot  be  carried  out  since 
the  layout  in  i  and  j  is  very  incomplete.  Suitable  explo¬ 
sions  have  taken  place  at  only  a  few  locations.  Asymmetries 
in  the  networks  and  the  scarcity  of  sources  means  that  (0<J-)ij, 
can  be  estimated  for  only  a  few  paths.  That  is,  complete 
travel  time  calibration  can  be  obtained  for  a  very  restrict¬ 
ed  set  of  source-station  pairs. 

However,  Model  1  has  a  more  serious  drawback.  The  tra- 
vel-time  tern,  T  <  ^  can  come  from  any  travel-time 

tables .  in  fact  one  could  use 

for  every  distance  and  depth  and  Model  1  would  still  apply. 

Thus,  in  this  respect  Model  1  is  trivial.  The  usefulness  of 
the  T(  Aij  hj)  term  is  that,  hopefully,  it  will  be  a  func¬ 
tion  of  only  the  distance  from  the  source  to  the  station 
and  will  be  on  the  order  of  a  few  seconds, 

now  consider  only  surface  events  and  set 


Let  us 
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T<A(/  ,0)  =  T(£uj). 

We  wish  to  define  a  .'world  wide  average  travel-time  curve" 
and  consider  to  what  extent  this  travel  time  is  a  function 
of  distance  only.  Computational  techniques  exist,  Engdahl 

et  al.,  (1968),  by  which  the  surface  focus  case  can  be  ex- 
tended  to  cover  nonzero  depths. 

Suppose  that  events  and  stations  can  be  placed  at  will 
throughout  the  world.  We  restrict  stations  to  be  on  surface 
land  masses  and  make  the  spatial  layout  in  i  and  j  complete 
with  the  requirement  that  all  distance  ranges  appear,  sup¬ 
pose,  further,  that  the  measurement  error  is  zero.  In  this 
case  the  model  for  an  observed  arrival  time  becomes: 

^m--o^T(Aa))  +  (0<s)-j({}  (5J 

where  i  refers  to  the  station, 
j  refers  to  the  source, 

(1)  is  a  dummy  subscript  and  refers  to  the  (fc)th 
distance  range, 

tij(|>)  is  the  observed  arrival  time, 

Oj  is  the  origin  time  of  the  source, 
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T(/^(Jl))  is  the  Predicted  travel  time  from  fixed 
tables , 

^  (t)  tlle  distance  from  source  j  to  station 

i,  and 

(0<S>ij(£)  is  an  interaction  time  term. 

The  interaction  term  accounts  for  the  failure  of  the  pre¬ 
dicted  travel  time  to  agree  with  the  observed  travel  time. 
The  model  (Model  2)  given  by  equation  (5)  describes  our 
population  of  travel  times.  That  is,  Model  2  gives  a  de¬ 
scription  of  any  specific  one  of  the  total  travel  times 
that  we  will  consider.  By  indefinitely  increasing  the  num¬ 
ber  of  sources  and  stations  the  population  approaches  the 
population  of  all  possible  travel  times  from  any  source  to 
any  station,  with  the  restriction  that  stations  are  on 
surface  land  masses. 

Known  geological  and  geophysical  properties  of  the 
Earth  indicate  that  different  regions,  will  have  different 
seismic  properties.  To  the  extent  possible  we,  further, 
assume  that  the  Earth  has  been  stratified  according  to 
differing  seismic  path  properties  and  that  an  equal  number 
of  events  occur  in  each  stratum.  We  stratify  primarily  to 
insure  that  the  population  gives  a  truly  "world  wide" 
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description  and  to  reduce  the  sampling  variance  of  any 

estimators.  That  stratification  can  significantly  reduce 

the  variance  in  sampling  is  well  known  (Cochran,  1963).  We 

consider  only  source-station  pairs  that  meet  the  distance 
requirement, 

Z5° s  £  too9 

and  suppose  that  each  event  is  seen  by  a  number  of  stations, 

V  WS  d°  not  re<l“ire  that  the  N.  are  all  equal.  We  will 
consider  this  further  when  we  discuss  randomisation.  We  will 
assume  that  all  events  are  in  the  magnitude  range  5%  to  7. 

We  make  this  assumption  in  order  to  keep  the  observational 
error  of  any  real  data  set  independent  of  the  randomisation 
errors  due  to  averaging  over  model  effects. 

With  these  restrictions  we  suppose  that  a  population 
Of  source-station  pairs  and  their  corresponding  travel  times 
is  selected  and  fixed.  This  population  of  seismic  events 
then  consitutes  the  population  of  events  for  which  we  will 
define  a  world  wide  average  travel  time  function. 

Further,  any  estimators  ultimately  considered  will  be 
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unbiasod  for  this  population  and,  in  general,  will  be 
biased  for  any  other  population.  This  emphasizes  the  care 
that  is  necessary  in  actually  selecting  the  population  or 
m  its  hypothetical  selection.  For  example,  in  estimating 
travel  time  corrections  (Tucker,  et  al.,  1968)  a  world  wide 
population  of  events  was  envisioned.  This  population  was 
suitably  stratified  and  events  selected  at  random  from  each 
stratum,  we  thus  believe  that  the  1968  Herrin  travel  times 
(Herrin,  et  al.,  1968)  are  unbiased  tor  this  world  wide 
population.  The  actual  problem  is  more  complex  than  we 
here  state  and  we  will  amplify  later.  However,  our  present 
point  is  that  the  sample  produces  unbiased  estimators  for 
the  population  sampled. 

Recalling  equation  (5),  we  see  that  each  travel  time 
in  the  population  can  be  expressed  as 

tym  +  (fiff)xj  K) . 

Define 


(6) 
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where 


and 


Aiu:joi  J  /N.. 

disianct  / 


00 


N..(H)  =  the  number  of  travel  paths  at  distance  %. 
Also  set 

V  _____ 

T  foff))  =  T(&(f))  +  (pa) " ,  (7) 

With  these  ■definitions  equation  (5)  can  be  rewritten  as 

ixjU)  0  °j  +  T*(At$+@<F%.Q' 


(8) 


(9) 


For  future  reference  observe  that 

0  foreacU  $ 

/JLti 

and  that  T*  is  a  function  of  only  the  distance 

from  the  source  to  the  event.  In  general,  (06)*  . 
varies  not  only  with  j  and  i  but  with  i  as  well.  However, 
in  view  of  <9)  this  is  a  local  change  with  &  for  each  j— i 
pair.  The  term  T*  (  ^  ^j)  is  the  world  wide  average  travel 
time  function  and  varies  only  with  |  (that  is  with  A  ). 

It  is  possible  to  give  a  physical  interpretation  to 
the  terms  in  equation  (8).  Suppose  that  the  Earth  is  a 
sphere  with  radius  r  (r  is  the  radius  of  a  sphere  with 
volume  equal  to  the  volume  of  the  Earth).  In  any  actual 
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computations  e  1 1  i  Dt  i  r^-S  4. 

P 1Clty  correcfcions  are  made  so  that  thi 
assumption  is  valid  „„ 

charact  • 

erased  by  a  velocity  distribution  V(f)  . 

velocity  distr-ihmt-  •  .  ’  P  given 

th,  „tl  M 

travel  time  t«  nce  the 

•  In  9®neral,  the  velocitv  .a 
ent  r  •  eiocity  structure  of  differ- 

-rons  of  the  Earth  differs. 

;7  aSSUmPUOnS  -  velocity  structure. 

Th°  re3i0"  C°ntaini^  -V  anomalous  velocity 

uctuae  as  confined  to  the  crust  and  upper  mantle. 

lateral  ^homogeneities  which  produce  these 

anomalies  are  confined  to  th 

°  the  uPPer  500-700  km  of 

e  Earth  and  below  that  depth  the 

pCi-  thG  velocity  dis¬ 
tribution  is  spherically  symmetric. 

(2)  The  velocity  distrihi.n 

trabutaon  of  any  re9ion  of  the 

Earth  is  such  that  for  distances  in 

Gs  m  excess  of  25° 

the  deepest  point  of  fv 

G  raY  path  is  below  500- 

700  km.  Thus,  the  p  =  i 

P  arrivals  that  make  up  the  pop- 

“  atl°n  Pa“  tWiCe  thr°^h  «-  upper  700  km  hut 
symmetrical  region  of  the 


mantle . 
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We  associate  T*  (  A  )  with  the  spherically  symmetrical 
portion  of  the  path.  Here  T*  ( A )  represents  a  virtual 
travel  time  dependent  on  distance  in  that  the  symmetrical 
portion  of  the  velocity  distribution  may  not  extend  to  the 
surface.  The  average  over  many  paths  represents  an  appar¬ 
ent  symmetrical  distribution  that  extends  to  the  surface 
and  which  would  produce  T*  (  A  ) .  Then  (0<S  represents 

deviations  from  a  complete  spherically  symmetrical  distri¬ 
bution.  Qualitatively,  we  see  that  tends  to 

decrease  with  increasing  A  since  the  inhomogenous  part 
makes  up  proportionally  less  of  the  total  path. 

Even  if  the  inhomogeneous  layer  extended  below  700  km, 
say,  included  the  entire  path,  the  interpretation  of  T*(A  ) 
as  a  virtual  travel  time  would  still  be  valid.  However,  in 
this  case  there  would  be  no  general  similarity  in  any  path 
m  the  Earth  and  the  usefullness  of  T*(A)  would  be  dimin¬ 
ished.  in  this  case  the  variation  in  (0G)*  would  be 
large  and  the  characterization  of  a  P  ray  path  in  terms  of 
epicentral  distance  would  have  little  meaning. 

We  will  henceforth  suppose  that  the  assumptions  (1) 
and  (2)  are  in  effect  and  that  T*  (A),  for  the  range  of 
A  of  interest,  is  of  the  order  of  hundreds  of  seconds  and 
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due  primarily  to  a  spherically  symmetrical  velocity  distri¬ 
bution.  we  assume  that  a  further  consequence  is  that 

l0<5>  ij(t>  13  at  most  an  °r!3er  of  magnitude  smaller  than 
T*  (A),  in  any  event,  T*  (A)  is  unique  once  the  population 
under  investigation  is  fixed.  Further,  for  an  increasing 
coverage  of  event-station  pairs  and  assumptions  (1)  and  (2) 

holding  the  dependancy  of  T*(A>  on  the  parent  population 
will  decrease. 

Equation  (8)  will  be  termed  Model  3  and  represents  the 
population  under  study.  Model  3  with  the  addition  of  obser¬ 
vation  error  will  be  called  model  4.  Model  4  is  not  trivial 
that  T*(£)  is  well  defined  and  essentially  unique.  Thus, 
Model  4  does  not  have  the  serious  defect  of  Model  1  but  again, 
in  practice,  neither  T*(A>  and  <0<J>*..(1)  can  be  estimated. 
Suitable  explosions  have  taken  place  at  only  a  few  locations 
and,  therefore,  the  layout  in  i  and  j  is  very  incomplete. 

We  also  define  Model  5  as  Model  4  without  the  depth  restric¬ 
tion.  Employing  the  computational  techniques  developed  by 
Engdahl,  et  al.,  (1968)  it  is  possible  to  determine  a  velo¬ 
city  distribution  from  T*(A).  Once  a  velocity  distribution 
is  established  it  is  possible  to  obtain  T*(A,h)  where  T*(A,h) 
represents  a  world  wide  average  travel  time  for  a  distance  A 


-14- 


and  with  an  event  at  depth  h.  It  is  important  to  note  that 
T*(^,h)  was  not  developed  by  "averaging"  as  was  T* (A ) . 
Thus,  T*(Zl/h)  can  only  be  interpreted  as  a  world  wide  aver¬ 
age  of  a  suitable  population  to  the  extent  that  a  spheri¬ 
cally  symmetrical  velocity  distribution  exists.  This  is  so 
since  the  computational  methods  of  Engdahl,  et  al.,  (1968) 
employ  T*  (A  )  and  its  derivative  to  "strip'1  the  crustal 
layers  from  the  surface-focus  curve  and  then  by  inversion  to 
obtain  a  velocity  distribution  which  is  employed  in  obtain¬ 
ing  T*(^,h)  for  h  nonzero. 

In  order  to  go  further  in  the  practical  analysis  of 
travel  times,  we  must  make  further  simplifying  assumptions. 
Since  assumption  (1)  seems  to  hold  in  practice  and  assump¬ 
tion  (2)  is  not  unreasonable,  if  we  consider  principally 
lower  order  effects,  we  shall  partition  the  interaction  term 
as  follows: 


(*a)ij(U  =  ftj(t)  ■-  GxjU) 


(10) 


where  is  primarily  associated  with  the  station  and 

6  .  .  is  primarily  associated  with  the  source, 

ij(D 


It  follows 


that 


(0<5)..a)  -  f/.'OD  + 


(in 
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where 


Am,  j  ai  J  /N 


dislanci  ( 


<*..({J=Z  6 a iu)/,  . 

MM  N-M) ' 

and  ls  as  before.  Then 

Wyg)  =  -(Jo),  _ 


(It) 


-  -f  G^j  -  _  Q" 


aj 


-  &j(ti  -  &.(£J  J  Gjjyy  _  Q  ^ 


m 


^%(t) +  Gij(fj 


(12) 


&jt)~  ?■■(£) 

a*jm=  G/-j(0~v..(J!) , 


and 
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Again  note  that 


Z  f **ja)  =  o 

All  i,  j  at  distance 


and 


Z  G'ijU)  =  ° 

All  i,  j  at  distance  Jt 


for  every  & 
=  Q  for  every  4. 


(13) 

(14) 


Our  models  now  become  Model  3a,  Model  4a  and  5a, 
respectively,  and  are  given  by 

=0j +_r  (ty))*  a)  + 

txjaik^ Qj  +T (  A» )4 $2^aJ  +  cr*.fl)  +  e fc  ,16, 


(15) 


and 


t 


°> fT  $/«J + a5y«)+ e ipk . 


(17) 


We  note  that  if  k  is  one,  that  is,  only  one  event  is  con¬ 
sidered,  then  it  is  impossible  to  separate  0* 

^ij(t)'°  i  j  ( Jt) 

ij(*)k*  ThUS  We  must  know  ^ij(|)  and  <5*iui)>  havo 
suitable  estimates  of  these  terms  or  have  a  simpler  model 

hold  if  we  are  to  observe  variances  of  the  order  of  0.01  to 
0.04  sec^. 

One  possible  simplification  was  made  by  Cleary  and  Hale 
(1966a),  who  assumed  that  the  partitioned  interaction  terms 
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^  i 3  (A)  and  ^  *i  j  (jl)  are  indePendent  of  distance  and  azimuth. 
The  resulting  "zero  th  order"  model  is, 


tjij  z  °i  +  T(^jj  hj h  +  6j  t  6/j 


(18) 


where 


0£  is  the  "station  time  term",  and 
ffj  is  the  "source  time  term". 

If  the  origin  time  is  unknown  then  Oj  and  6.  are  statistically 
confounded  so  that  we  may  write 


t 


*j  =  0j  +  T*(A,j>j)  +  ^  +  e-. 


(19) 


where  0*  is  the  "biased"  origin  time. 


It  is  possible  to  express  the  zero  th  order  model  terms 
0i  and  ffj  as  suitable  means  of  and  <5*ij(Jl)  and  the 

error  term  6  as  a  sum  of  and  mean  deviations  in 

the  andQ*ij(|)  terms  of  Model  5a.  However,  it  is 

impossible  to  give  a  randomization  interpretation  if  k  is 
one  so  that  we  will  not  pursue  this  further.  We,  thus,  inter¬ 
pret  and  CTj  as  true  model  terms  that  are  not  functions  of 
distance.  Furthermore,  the  error  term  €i;j  is  also  assumed 
to  be  independent  of  distance.  In  fact  £  L  .  is  just  measure¬ 
ment  error  under  the  models  given  by  equations  (18)  and  (19). 
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If  the  zero  th  order  model  is  adequate,  the  sum-square 
of  residuals  (suitably  normed)  resulting  from  locations  of 
earthquakes  used  to  estimate  0j*,  T*<  A,h)  and  ^  should  be 
comparable  to  the  variance  of  £...  In  the  cleary-Hales  (1966) 
studies  this  sum  was  of  the  order  of  1  sec2  (Hales,  personal 
communication,  1969),  much  larger  than  expected  from  the 
simplified  travel  time  model  (less  than  0.1  sec2),  clearly 
ne  error  term,  €iy  in  equation  (19)  was  absorbing  more 
"error"  than  would  be  expected  due  to  uncertainties  in  measure¬ 
ment  and  timing  alone.  The  Cleary-Hales  studies  indicate  that 

the  simplified  model  is  inadequate  in  explaining  the  interaction 
terms  of  Models  3a,  4a  and  5a. 

Later  Douglas  and  Lilwall  (1968)  made  a  further  model 
simplification  in  studying  essentially  four  events  in  the  Pacific 
whose  locations  were  known  from  independent  studies.  Douglas 
and  Lilwall  in  effect  assume  that  <5  ^  is  constant  for  the  four 
event  areas.  The  data  are  from  the  Eniwetok  series,  Bikini  series, 
a  large  Hawaiian  earthquake  and  Longshot.  The  standard  errors 
of  observation  as  quoted  by  Douglas  and  Lilwall  are  Longshot, 

0.8  sec,  Hawaii,  1.2  sec,  Eniwetok:  0.  5.  sec  (average)  and 

0.5  sec  (average).  The  range  on  Eniwetok  is  0.3-0. 7  sec 
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and  on  Bikini  0.4-0. 6  sec.  Thus  it  again  appears  that  the 
error  term  in  a  simplified  model  is  absorbing  more  "error" 
than  would  be  expected.  In  this  model  the  station  time  tex'ms, 
0^,  are  independent  of  source  location  and  could  thus  be 
equally  well  determined  (i.e.  unbiased  estimates  obtained) 
from  any  given  explosion.  Joint  estimation  using  four  events 
serves  only  to  reduce  the  variance  of  the  estimates  under  the 
assumed  model. 

It  follows  as  an  implication  of  this  model  that  estimates 
of  station  time  terms  obtained  in  the  manner  just  described 
should  be  equally  applicable  to  sources  anywhere.  Ignoring  a 
base  line  shift  (which  could  be  due  to  not  estimating  the  con¬ 
stant  source  term)  the  1968  Douglas-Lilwall  estimates  differ 
from  those  given  by  Lilwall  and  Douglas  (1970).  A  further 
consequence  of  this  model  is  that  the  estimated  time  terms 
should  never  appear  to  be  a  function  of  source  location.  That 
is,  time  terms  must  not  appear  to  vary  with  distance  or  azimuth 
from  station  to  epicenter,  Cleary  and  Hales  (1966b)  later  re¬ 
ported  that  in  fact  "the  remaining  residuals  were  azimuthaliy 
dependent  for  some  stations".  Also  Lilwall  and  Douglas  (1970) 
obtained  azimuthally  dependent  estimates  which  for  many  stations 
are  statistically  significant.  Thus,  the  simplified  model 
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and  on  Bikini  0.4  0.6  sec.  Thus  it  again  appears  that  the 
term  in  a  sirop.lified  model  is  abr  orbing  more  "error” 
than  would  be  expected,  in  this  model  the  station  time  terms, 
0.i»  are  independent  of  source  location  and  could  thus  be 
equally  well  determined  {i.e.  unbiased  estimates  obtained) 
from  any  given  explosion.  Joint  estimation  using  four  events 
serves  only  to  reduce  the  variance  of  the  estimates  under  the 
assumed  model. 

It  follows  as  an  implication  of  this  model  that  estimates 
of  station  time  terms  obtained  in  the  manner  just  described 
should  be  equally  applicable  to  sources  anywhere.  Ignoring  a 
base  line  shift  (which  could  be  due  to  not  estimating  the  con¬ 
stant  source  term)  the  1968  Douglas-Lilwall  estimates  differ 
from  those  given  by  Lilwall  and  Douglas  (1970).  a  further 
consequence  of  this  model  is  that  the  estimated  time  terms 
should  never  appear  to  be  a  function  of  source  location.  That 
is,  time  terms  must  not  appear  to  vary  with  distance  or  azimuth 
from  station  to  epicenter.  Cleary  and  Hales  (1966b)  later  re¬ 
ported  that  in  fact  "the  remaining  residuals  were  azimuthally 
dependent  for  some  stations".  Also  Lilwall  and  Douglas  (1970) 
obtained  azimuthally  dependent  estimates  which  for  many  stations 
are  statistically  significant.  Thus,  the  simplified  model 
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appears  to  be  inadequate  in  that  the  interaction  terms, 

^ij(l)  and  <3*ij(|)'  can  not  adequately  be  represented  by 
two  simple  time  terms. 

Cleary  and  Hales  (196&)  and  Herrin  and  Taggart  (1968a) 
have,  using  statistical  analyses  of  earthquake  data,  shown 
that  the  station  time  terms  for  a  number  of  stations  are 
dependent  upon  azimuth  from  station  to  source.  These  find¬ 
ings,  however,  may  not  be,  strictly  speaking,  pertinent  to 
this  discussion  in  that  they  are  based  on  a  simplified  travel 
time  model.  To  use  these  results  as  evidence  for  keeping  the 
interaction  term  in  the  travel  time  model  is,  in  a  sense, 
begging  the  question.  However,  there  are  other  studies  of 
azimuthal  effects  which  do  not  suffer  from  this  drawback. 
Ritsema  (1959)  showed  that  normalized  amplitudes  of  p-waves 
recorded  at  Djakarta  varied  sinusoidally  with  azimuth.  Using 
P  and  p?  phases  recorded  by  the  Berkely  network,  Otsuka  (1966) 
found  the  time  terms  to  be  sinusoidal  functions  of  azimuth  to 
the  event.  Nizai  (1966)  reported  the  same  effect  at  the 
Tonto  Forest  Array.  Bolt  and  Nuttli  (1966)  studying  data 
from  the  Berkely  network  showed  conclusively  that  the  time 
terms  for  Shasta  and  Marysville  varied  sinusoidally  with 
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a^uth.  The  amplitude  of  thc  variation  wss  grMter  ^  i 
Errors  in  epicentral  determinat  ^  ^  ^  ^ 

time  tables  used  enter  onlv 

nly  as  second  order  effects  in  these 

studies. 

Reciprocity  retires  that  source  as  well  as  station  ^ 

MSt  ^  C°nSidered'  in  to  vary  with  asimuth 

-  an  explosion  ware  fired  beneath  the  Shasta  station  we  „ould 

expect  to  see  a  variation  of  travei  time  with  asimuth  of  at 

least  ±  l  sec  at  teleseismic  distances  Thus 

'  Thus  we  are  forced 

return  to  the  more  complicated  tr*  vpi  • 

r  xuuuea  travel  time  model, 

With  mild  regularity  assumptions  we  can  expand  the  two 

coefficients  being  functions  of  epicentral  distance.  Thus, 

M*  -A  ^  •  - 


•  i jit)  *  I  fyf  S,kl  +  i  ?• ' 

J 

+  2  sin  22y.  +  2  B u coS t 

J 

s*^d)  <5<  +  ,  flj/ s'"  ty  +  ,  By*  c« 


+2  Qjjt  sih  zfy  +  z  By* 
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where  "1"  indicates  the  dependence  of  a  parameter  on  epicentre  1 
distance,  AU)  ■ 

We  now  write  a  first  order  model,  neglecting  second  and 
higher  order  terms  in  the  Fourier  expansion,  as  follows: 

+  T*^-)  i  <pxi  i 
+  <osZ^  +  e*^  (20) 

where  ztj  is  the  azimuth  angle  from  station  to  source,  z 
is  the  angle  from  source  to  station  and  6*.  is  e  .  . 

i](l)  ls  C13(1) 

plus  the  remaining  Fourier  terms  of  the  expansion  of  0* 

and  lj(I) 

C  •  Note  that  the  symmetry  of  doubly  indexed 

variables  has  been  lost  in  this  expression,  not  because  of 
a  failure  of  reciprocity,  but  because  the  two  azimuth  angles 
associated  with  the  ends  of  a  great  circle  arc  are  not  simply 
related.  This  travel  time  model  is  an  expansion  to  the  first 
order  in  azimuth  angle,  of  the  complete  interaction  model 
(Model  5)  under  the  assumption  that  the  interaction  term  can 
be  separated  into  a  source  effect  and  a  station  effect. 

Using  the  explosion  data  where  0j/  A (1)  and  hj  are 
known,  we  can  attempt  to  estimate  the  parameters  in  the  first 
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order  model,  provided  sources  and  receivers  are  well  distri¬ 
buted  geographically  and  in  distance  ard  azimuth  relative  to 
one  another.  .Attempts  to  make  such  estimates  have  generally 
failed  because  known  sites  of  large  explosions  are  rare. 

A  further  simplification  can  be  obtained  if  the  dis¬ 


tance  dependence  of  the  Fourier  coefficients  is  neglected. 
As  final  model  forms  we  consider,  Model  6, 

"  Pj  +  4  *4  ft l  Stn  B  ' ' 

sj 

+  B -Kq-KijsIfiZji 

i  bj  c osZji  +  £**(q  on 


and,  Model  7, 


Uj(tf  Oj +  $1 +  A; 

tB: 


-j  +  ajsinty 
+ b< COJ  + e!$J 


where  +  is  plus  the  difference  between 

and  0i+  AiSinZij  +  B-jcosZ^  and  the  difference  between 

®  *ij  (1)  anfl  Si  +  aj  sin  zji  +  »>j=os  Zji  . 

Models  6  and  7  represent  expansions  to  the  first  order  in 


azimuth  angle  with  constant  coefficients  under  the  assumption 
that  the  interaction  term  can  be  separated  into  a  source  effect 
and  a  station  effect. 


STATISTICAL  PROBLEMS  IN  THE  ESTIMATION 


OF  MODEL  PARAMETERS 


With  the  exception  of  measurement  errors,  the  popula¬ 
tion  described  by  Model  5a  is  deterministic;  that  is,  repeated 
observations  on  a  fixed  travel-time  path  will  result  in 
essentially  the  same  travel  time.  Any  simplified  model,  for 
example,  Model  7,  will  thus  have  a  fixed  error  term  that 
represents  the  difference  in  the  true  travel  time  and  that 
predicted  from  the  simplified  model.  In  an  attempt  to  treat 
this  fixed  error  term  we  assume  that  any  seismic  event  can  be 
considered  to  have  been  selected  at  random  from  the  population 
of  possible  events.  Observe  that  if  we  do  not  employ  random¬ 
ization  that  this  model  difference  remains  a  fixed  and  unknown 
time  error.  Hence  any  estimation  techniques  must  be  arbitrary. 
To  avoid  this  dilemma  we  chose  to  employ  randomization. 

Consistent  with  our  model  studies  we  suppose  that  a  fixed 
set  of  events  and  stations  has  been  determined.  How  repre¬ 
sentative  of  the  whole  Earth  our  final  results  are  depends  on 
the  population  of  source-station  pairs.  Thus  the  1968  P 
travel  times  may  be  biased  toward  the  seismically  active 
regions  of  the  Earth  since  approximately  90%  of  the  data 
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employed  came  from  these  regions.  Also  the  results  may  be 
biased  toward  earthquakes  since  only  13  of  the  total  of  292 
events  employed  were  shots.  That  is,  we  cannot  reasonably 
assert  that  the  sample  employed  in  the  1968  P  travel  time 
study  (Special  Number  -  1968  Seismological  Tables  for  p 
Phases,  Bull,  seism,  Soc.  Am.  58  No.  4)  was  a  random  sample 
from  a  world  wide  population  or  a  population  with  extensive 
shot  data.  However,  if  we  restrict  the  scope  of  the  popula- 
ion  then  we  can  claim  that  we  obtained  a  random  sample,  it 
is  to  this  restricted  population  that  we  make  inferences. 

Until  data  are  obtained  in  the  less  active  regions  of  the 
Earth  there  is  no  way  to  assertain  the  usefulness  of  the  1968 
P  travel  times  for  use  in  these  areas.  However,  in  view  of 
the  known  geological  and  geophysical  properties  of  the  Earth, 
we  believe  that  the  1968  times  should  be  quite  useful;  recall 
that  a  spherically  symmetric  velocity  distribution  appears  to 
hold  for  the  Earth.  Thus  we  believe  that  the  restricted  popu¬ 
lation  is  representative  of  the  whole.  Earth  as  regards  earthquakes. 

We  could  similarly  argue  that  the  results  would  apply  to 
shots  in  general.  In  view  of  independent  work  by  Veith  and 
Clawson  (1972)  we  shall  consider  the  shot  data  problem  somewhat 
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differently.  Figure  1  and  Tabl  , 

nl  1  are  excerpted  fret.  Veith  and 

Clawson  (1972)  which  the  authors  were  hind  „  „ 

.  „  fcu,d  enouSh  to  let  us  use 

before  publicatior  p-:~,  , 

*  Flgure  1  contains  a  plot  of  fha 

pj-ot  of  the  residuals 

°”  1968  P  tiBeS  °f  «“  given  in  Table  1  In 

preparing  figure  1  multiple  ohservations  of  travel  times  along 

the  same  path  were  averaged  and  considered  as  one.  also  any 

before  plotting  in  Fiqure  1  mv 

s  figure  ±.  The  removal  of 

01  ^ne  mean  source 

effect  will  reduce  the  variability  in  the  data  but  should  not 

13  18  the  °“  ~  communication.  1971).  Furthers 

no  statron  corrections  were  applied.  Thus  Figure  1  presents 
-  residuals  fr..„  the  19SB  p  times  in  .  ^  „f  43  ^  ^ 

-44  observations  at  teleseismic  distances.  The  data  set  should 

not  be  biased  due  to  path  overloading  or  mean  removal.  The  re- 

suits  are  most  pleasinrr  t,.:  . 

sual  inspection  reveals  that  there 

13  ^  aPPreCiablS  ““  the  zero  axis  beyond  25. 

.  In  °rder  t0  quantitatively  investigate  the  data  the  actual 
residuals  were  read  from  Figure  1  by  overlaying  graph  paper. 

The  readings  were  made  to  .1  second  mi, 

second.  The  results  of 

anyalyzing  this  data  on  a  5"  cell  h„  • 

cell  basis  are  presented  in 
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Table  2.  Due  to  the  close  agreement  of  the  computed  standard 
deviations  with  the  plotted  curve  calculated  on  a  2°  basis 
we  believe  not  much  accuracy  was  lost.  Again  the  results  are 
pleasing.  In  the  distance  range  25°  -  100°  the  largest  mean 
is  .192  sec.  at  25°  -  40°.  The  largest  relative  error  is 
approximately  one  part  per  thousand  and  occurs  in  the  30°  -  35° 
and  35°  -  40°  cells.  These  errors  are  so  slight  that  they 
would  never  affect  a  seismic  location  with  reasonable  station 
coverage.  In  fact  if  one  is  only  concerned  with  location  then 
the  removal  of  the  source  mean  from  the  residual  data  has  no 
effect  since  for  adequate  networks  this  could  only  bias  the 

origin  time  estimate  (Tucker  et  al.,  1968). 

If  we  assume  that  the  observations  are  independent  normally 

distributed,  the  efficacy  of  which  can  be  questioned,  then  we 
can  make  t  tests  for  each  cell.  The  t  values  are  in  the  last 
column  of  Table  2.  If  we  choose  an  °<  of  .01,  then  the  means 
of  cells  25°-30°,  30°-35°  and  95°-100°  are  significantly  dif 
ferent  from  zero.  The  degrees  of  freedom  for  any  of  these  te„ ts 
are  so  large  that  they  are  effectively  infinite.  The  critical 
t  value  is  2.576  for  a  two-sided  test.  Thus  only  the  35°-40°  . 

cell  appears  to  be  definitely  significant.  However,  since  any 
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travel-time  curve  must  possess  certain  differential  properties 
we  believe  a  better  procedure  is  to  fit  a  smooth  curve  to 
the  data.  With  randomization  we  can  obtain  -optimal"  estimates 
of  the  fitted  curve  without  the  assumptions  of  independence 
and  normality,  we  return  to  these  considerations  after  we 
finish  the  discussion  of  the  randomization  procedure. 

We  thus  believe  that  the  1968  P  times  are  not  seriously 
biased  against  the  shot  data  population  of  the  Vei th-clawson 
work.  Again  geological  and  geophysical  considerations  lead 
us  to  believe  that  the  restricted  population  is  representative 
of  the  whole  Earth.  We  conclude  that  while  we  sampled  a 
restricted  population  in  the  1968.  studies  that  this  population 
is  representative  of  the  whole  Earth  with  respect  to  both 
earthquakes  and  shots,  we  will  consider  this  point  again  in 
the  section  on  estimation. 

However,  representative  or  not,  we  envision  a  final 
population  of  event-station  pairs  and  their  corresponding 
travel  times,  we  list  these  events  as  shown  in  Figure  2  and 
suppose  that  there  are  K  strata,  M  events  and  N  stations  in 
the  population.  Hot  all  stations  see  all  events,  that  is, 
are  within  the  distance  restrictions.  Also  each  event-station 
pair  is  at  some  unique  distance  value. 


We  suppose  that  the 
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range  of  A  has  been  broken  into  cells  and  that  there  are  L 
of  these.  We  assume  that  the  cells  are  of  equal  width  and 
small  enough  that  the  effects  of  grouping  are  negligible. 

The  events  are  classified  first  according  to  stratum  and  second 
according  to  event  within  a  stratum,  that  is,  all  events  are 
different  so  that  event  1  in  stratum  1  is  different  from 
event  1  in  any  other  stratum  and  so  on.  However, the  N  stations 
under  each  individual  event  are  the  same  N  stations,  that  is, 
station  1  in  stratum  1  -  event  1  is  the  same  station  1  in  all 
other  stratum-event  combinations.  since  each  event-station 
pair  has  a  unique  distance  there  is  one  travel  time  value  for 
each  row  of  the  figure.  These  are  indicated  by  the  x's  in  the 
body  of  the  figure.  The  randomization  scheme  is  given  by  the 
following  procedure: 

(1)  Select  by  simple  random  sampling,  that  is,  with 
equal  probabilities  and  without  replacement,  k^  of  the  K 
strata  where  li  k^  ±  K. 

(2)  Select  by  simple  random  sampling  n^of  the  M  events 
in  each  stratum  selected  in  (1)  where  l*jm^«,M. 

(3)  Let  be  the  number  of  stations  that  see  the  mth 
event  in  the  kth  stratum.  Select  by  simple  random  sampling 

of  the  Nj,m  stations  for  each  event  selected  in  (2)  where 
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1  a  "km  *  Nk >  • 


the  1968  studies  steps  (1)  and  ,,, 

’  and  (3)  were  carried  out 
(Tucker,  et  al  iqcq  ,  Ut 

•'  1968  an<*  Herrin,  et  al  iq.o, 

data  were  1968) ‘  The  available 

were  grouped  on  a  1°  latitude  bv  1»  , 

the  laro  t  longitude  basis  and 

he  largest  source  event  selected  e 

a  i  fr°m  EaCh  ^oup.  The  groups 

are  the  strata  and  k,  =  K  „  P 

1  K.  We  assume  that  tv.o 

i  tnat  the  largest  event 

not  behave  systematically  in  the 

so  that  tv  7  he  Population 

the  se  let  ion  of  events  i„  the  strata  c  v, 

strata  can  be  thought 

lent  to”"  °”‘  -Here  M  =  1  and  MCh  COmbinati°n  iS  •*>*.- 

J  to  a  particuiar  ,  in  our  previous  notation.  SteP(3)was 

1968  that  nk”  =  Kk”’  The  Seri°US  dlfficulty  with  the 

1968  procedure  is  that  the  cell  or 

with  equal  •  Pln9  WlU  n0t  Pr°dUce  stnata 

shVZ  ;elsmic  re91ons  the  assumption  of  —  -  — 

should  be  adequate.  Note  thet  • 

iMuce  that  since  n,  —  »T  ,  n 

Nkm  we  have  km/Nv  =1 
Now  consider  hm  X’ 

onsiaer  the  general  randomization 

luomization  scheme  with  the 

restriction  that  nkm/Nv  =  Of  , 

'km  /  a  constant  not 

.  '  necessarily  equal 

-  -  ™e  can  show  that  the  probability  that  an  ind  vidua  1 

event-station  pair  is  included  in  a  sample  is  (Wl/M)o« 

Tehran,  1963K  Thus  the  inclusion  probabilities  are  equa, 

consider  a  given  distance  ranqe  th  S 

9  •  hen  the  number  of  observations  in 
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the  whole  sample  at  this  range  is  a  random  variable,  sim¬ 
ilarly  the  number  of  observations  made,  by  any  station  is 
random.  Thus  it  is  not  obvious  that  the  usual  statistics 
(sample  mean,  sample  variance,  regression  estimates,  etc.) 
are  applicable. 

We  shall  now  discuss  a  general  estimation  scheme  that 
will  justify  the  use  of  the  ordinary  statistics  and  indicate 
some  further  implications  of  randomization.  The  technique 
is  a  generalization  of  the  linear  regression  estimator  in 
sampling  theory  (Cochrarj  1963).  Suppose  we  have  a  population 

of  observables  Yi#  i  =  1,  2 . N.  We  would  like  to 

predict  the  values  in  terms  of  a  set  of  independent  vari¬ 
ables  Xxi' x2i'  •••'  >  that  is,  for  given  i 

and  hence  Xli#  ...,  Xj^  we  desire  an  estimate  of  Y^.  In 
order  to  obtain  a  solution  we  assume  that  the  Y’s  can  be 
suitably  estimated  by  linear  combinations  of  known  functions 
of  the  X's,  that  is,  we  assume  a  set  of  basis  functions  exists. 
We  employ  these  basis  functions  in  only  finite  sums  so  that 
we,  in  genera],  cannot  approximate  the  Y's  exactly.  We  write 


=f><>  4 Pi hftuj 


23) 


we  suppose  that  the  |j'j  are  such  that 

M  , 

QC(3 )  =  L  Sj 

I  ^~i  A 

is  a  minimum.  Thus  the  (population)  predictor  is  determined 
by  least  squares. 

This  definition  of  the  predictor  is  admittedly  arbitrary 
but  can  be  justified  as  follows,  suppose  that  an  individual 
Ya  is  selected  at  random  from  the  population.  Then  one  by 
randomisation  theory  (Kempthorne,  1955;  cornfield,  !944)  can 
show  that  if  we  express  this  y.  as  in  equation  (23)  with 


)  (24) 
mean  and  can  be 


then  becomes  a  random  variable  with 
expressed  as 

-  H  +  e;  • 

The  mean  of  is  zero  and  the  variance  of  d 

A /  . 


(25) 


is 
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we  see  from  equation  (26)  that  V(  «;)  is  a  minimum  consistent 
with  our  model  requirements,  that  is,  linear  basis.  If  no 

Yi  is  any  more  important  than  any  other  than  our  assumptions 
on  then's  is  not  unreasonable. 

Since  the  population  predictor  is  based  on  linear  least 

squares  we  have  an  immediate  solution  for  then’s  (Graybill, 
1961)  ; 


^(x,x)~'x'yj 


where 

I  W  -A.)  ••• 

I  ^1  >  •  - 
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and  denotes  inverse  and  *  a 

fl  's  are  tra"Spose-  ™us  the 

P  -  ilnear  functions  Qf  the  pQpu 

it  is  often  imr*  -V,  In  practice 

^n  impossible  t<-> 

to  determine  the  fl  •=  „ 

know  the  y.  ,k  P  Slnce  we  0°  not 

c  x  s.  Thus  we  mnc4-  • 

must  estimate  the  A  . ,,  . 

e.  P  s  by  some  scheme 

^  (i's  are  linear  in  the  y's  if 

random  sample,  then  „  0  ^  '  *  "*  take  a  simpie 

optimum  estimate  of  X'Y  i<=  • 

W..  M../,.  X  V  is  given  by 


where 


=fSampfet|  rows  of  )(J 

^m=  (SamplaJ  rows  of  Vl 

^  -  is  the  sample  site,  that  is,  x  is  th 

basis  functions  and  y  is  the  „  SamPled 

The  esti  t  a"  C°rreSPOndin9  «*«•  V  values, 

me  eotimator  X'Y  it;  mu  j 

mat  6  3  Horvit 2 "Thompson  (ht)  esti¬ 
mator  and  is  onlv  *.  •  esti- 

unbia  d  ■  eStimat°r  SlnCe  "0  uniformly 

-ed  minimum  variance  estimator  of  the  totals  „  exists 

sampling  from  finite  popuiations  (Godambe,  1955'  fjan 

1968)  yoo,  Hanurav, 
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t  complete  properties  see  « 

see  Hanurav  (1968). 
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Thus 


E(x'y)*x'y 


so  that 


E 


P 


as  we  see  from  equation  r?7\  mi 

(27).  Thus  an  ont  i  mn  i  .  • 

P  optimal  estimator  of 

is  given  by 


Hx'x)>)y'yn . 


(28) 


In  order  to  employ 
In  the  1968  studies  thes 
we  have 


equation 


(28)  we  must  know  n  and  x 


®  quantities  were  not  known. 


However, 


frw'feO-  (S  x'xj1 
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,  .  _  v,  /V  y  . ) .  Tf  the  sample  chosen  is 

where  hj(x)  =  hj(xii'***'  *kiJ.  a*  ^ 

representative  of  the  population  X's,  then 


X*X  =  Sx'x 


(29) 


in  fact,  if  the  strata  are  such  that  the  X's  are  grouped  into 
m  equal  size  sets  of  common  values  and  we  sample  such  that 
one  from  each  common  set  is  chosen,  then 

XmX«>=  n'X'X  i  x 

as  one  can  easily  show,  unfortunately  this  was  only  approx 
imately  done  in  the  1968  studies  since  not  all  events  see  all 
distance  ranges  and  there  was  no  stratification  with  respect 

to  stations. 

However  replacing  (m/N  X'X)'1  by  (X'm  X,*)'1  is  analogous 
to  forming  a  linear  regression  estimator  (Cochran, 1963) .  We 

thus  consider 

^=(x:xJx:y 

which  is  the  usual  linear  least  squares  estimator  (Graybill, 
1961) . 

We  will  now  develop  a  bias  formula  for  the  general 
estimator  given  by  equation  (30).  From  equation  (25)  we 


y„=x„p 


m 


have 


(31) 
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where 


-  (sa„plj  e's  J 


Thus 


-  P  +(XmV^,)  'Xto  e„ 


v\  , 


e(x'e.)=SX'c 


where 


e./i 


But  from  equation  (25) 


(32) 


x'c'X'(Y-Xt) 

=*'y-(x'x)P 

= y'y  -  x'y 

=  <t> 

Thus,  in  repeated  sampling  with  fi*  a  •  <34) 

9  lth  flxed  size  m,x'  p  4.  , 

random  vector  distributed  about  a  zero  m 
,  Zero  mean  vector  and  for 

“ge  m  the  »i«  should  be  small. 

We  have  from  equation  (32)  that 
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if  rLyV-^i  XX  in  probability  and  m  is  suitably 

f' Ww'm  fsj 

large.  From  equation  (36)  we  obtain,  neglecting  higher 
order  terms, 


mu  /  | 

(f-  p>  i  ft  *'*)"{;  O^x'xyl'e. . 


(37) 


Let  the  entries  in  any  X  matrix  be  denoted  by  x. .  and 

•**  J 

, 

the  entries  in  (X'X/N)  1  be  denoted  by  X^j.  Set  p=n  +  1. 


Then 


x/^ftx'x)‘x,X 


5  US  t 


Vv\ 


XJ 
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aS  one  maY  5how  bY  straight  forward  but  tedious  algebra. 
We  have  for  the  qth  entry  in  equation  (38)  that 
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where 

and 


f  =  m/N 


ii  «v, 

S„=  Z  Z  2  Xjr.  "Kb*: 

°  j£i  j=i  J  n  A 


pX'^'  /  V 

J  4n-2  * 


(40) 


(41) 


Finally,  from  equation  (37)  we  have  that 

Therefore  the  bias  is  of  order  1/m  if  Sk  is  of  order 
unity  for  every  k.  A  sufficient  condition  for  this  is  that 
the  and  Yn- i  c  be  bounded.  This  is  the  case  for  the 

X  J  S  x  o 

data  of  the  1968  studies.  Thus  the  bias  due  to  randomiza¬ 
tion  is  of  order  1/m  for  these  studies. 

Following  Cochran  (1963  pages  194-1^)  one  can  show  that 

an  a s symptot ica lly  valid  estimator  of  , 2-  S’ 

l  * 

{ltd'll  [I  X/ 

4m- k) 


13 


(42) 
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the  usual  estimator.  Recalling  equation  (30)  we  see  that 

the  usual  least  squares  estimators  are  appropriate  with  our 
randomization  model. 

Thus  employing  the  usual  least  squares  estimators  should 
produce  negligible  bias  in  the  travel  time  estimates  since  in 
the  distance  range  25“  -  100°  the  fewest  observations  were 
81  in  the  100°  cell.  The  case  for  station  terms  is  not  as 
conclusive  since  we  reported  estimates  for  stations  with  as 

few  as  10  observations.  However,  most  of  the  station  terms 
were  estimated  with  much  more  data. 

While  we  have  assumed  simple  random  sampling  in  proving 
these  results,  the  HT  estimator  is  optimum  under  general 
randomization.  Also,  if  the  design  is  balanced  so  that  the 
inclusion  probabilities  are  equal  and  further  the  second 
order  inclusion  probabilities  are  equal,  then  the  bias  is 
again  of  order  1/m.  We  have  previously  shown  that  the 
(first  order)  inclusion  probabilities  are  equal.  Since  we 
sample  all  strata  and  use  all  stations  that  see  any  event 
the  second  order  inclusion  probabilities  are  either 
1,  m (m-1 ) /M (M-l )  or  (m/M)2 .  We  obtain  1  when  we  consider 
two  stations  that  see  the  same  event,  m (m-l )/M (M-l )  when 
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WS  C°nSlder  two  events  in  the  same  stratum  and  (m/M) 2 
when  we  consider  two  events  in  different  strata.  Thus 
the  second  order  inclusion  probabilities  are  not  equal. 
However,  from  the  derivation  of  equation  (39)  it  is 
Obvious  that  with  these  inclusion  probabilities  again 
the  bias  is  of  order  1/m.  Mso  with  a  sufficiently  large 
total  sample  size  the  usual  variance  estimator  is  also 
applicable.  Thus  these  results  apply  i„  the  randomization 
scheme  employed  in  the  1968  studies  and  the  usual  esti- 
mators  are  applicable  in  that  work. 

bet  us  now  consider  the  general  case  in  which  measure¬ 
ment  error  is  also  present.  Thus  suppose  that 


where 


Kv\  V"  p  +  ^  ry^  +  £ 


e  ~N($j  <52I) 


(43) 


Then 


whrch  rs  the  result  given  earlier  in  equation  (32). 
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Also 


f  (y»  [l-XJtixj'x*]  X»/  , ) 

4m  -  k) ' 

=(  A)  [ E  ^'li-x^xx 

+  E  e/[l-)C^Xd~lX'mJe? 

since  6  and  are  independent_  ^  ^  ^ 

we  have  that 

E  (X-p+eJtl-XUx^Xj'xJ: ](Xm)/  , 

.  w/  ,  /fm-k) 

-  Vte). 

Also  it  is  well  known  (GraybiU,  1963)  that 

Ee  [  I'XVXm  Xon)  Xfie/Ctn_Q 


Thus 


=  <s2. 


4x » [i  -  x.r  'x:j  y,/^  t) ) 

=  V6*u)  +  c2 

and  again  the  usual  estimator  is  appropriate.  Finally, 
in  the  general  case  the  usual  least  sguares  estimators  are 
symptotically  unbiased.  If  we  require  that  the 
regression  estimator  be  optimal  with  .respect  to  the 

randomisation  scheme  then  the  derived  estimator  is  the  only 
one  since  the  HI  estimator  is  unique.  ' 
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We  now  consider  one  final  result  of  randomization. 
Suppose  that  a  single  event  is  drawn  at  random;  that  is, 
one  stratum  is  selected  at  random  and  then  one  event  is 
selected  at  random  within  the  selected  stratum.  We 
suppose  that  Model  4  holds  and  is  known.  With 

randomization  becomes  random.  Also  suppose 

that  for  each  event  in  the  population  each  station  is  at 
a  different  distance  cell.  Then  for  a  fixed  distance  5 

l  ^  > 

we  have 

£(  E  HjYiJ  +  e) 


=  «) 
fiUyH  J 


=  o 


from  equation  (9) .  Thus,  if  we  draw  random  samples  from 
our  population,  we  expect  that  at  each  distance,  1 ,  the 
residuals  from  t|A)  would  average  out  to  zero.  This  is 
essentially  the  case  with  the  Veith-Clawson  data.  It  is 
in  this  sense  that  a  world-wide  average  travel-time 
function  is  unambiguously  defined. 


Note  further  that  if  an  event  is  fixed  then  the 
i  j  ( )  are  not  random.  The  only  random  quantity  is  {£ 
If  one  is  then  attempting  to  locate  this  fixed  event  with 
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a  fixed  set  of  stations,  then  the  location  may  be  signifi¬ 
cantly  biased.  Longshot  is  an  example  of  this  phenomenon. 
If  repeated  shots  were  made  on  Amchitka  and  a  given  set  of 
stations  used  to  locate  these  shots,  then  the  1968  travel 
times  would,  in  general,  give  locations  approximately  25  km 
to  the  north.  One  could  argue  that  the  1968  times  are  not 
"tuned"  to  Amchitka  and  this  is  the  case.  However,  if  the 
1968  times  are  altered  to  fit  Amchitka,  the  altered  times 
would  not  fit  the  Veith-Clawson  data.  We  shall  return  to 
this  point  when  we  consider  possible  models  of  the  Veith- 
Clawson  data.  This  being  the  case,  it  should  be  obvious 
that  how  well  a  given  world-wide  average  travel-time 
function  locates  a  specific  event  is  not  a  suitable 
criterion  of  the  validity  of  the  travel-time  function. 

The  problem  now  arises  of  how  to  make  statistical 
tests  with  the  randomization  scheme  employed.  Three 
approaches  are  possible.  Oris  is  to  consider  non-parametric 
techniques.  We  rule  these  out  since  ranks,  runs,  etc.  of 
travel  times  make  little  sense.  A  second  possibility  is 
to  employ  as  symptotic  theory  to  justify  the  normality 
assumption.  In  general,  this  requires  conditions  on  the 
higher  order  moments  of  the  population  (Madow,  1948). 


In  our  problem  it  is  impossible  to  verify  these  conditions. 
However,  with  "reasonable"  sample  sizes  for  many  problems 
the  assumption  of  normality  is  adequate  (Kempthorne,  1955) . 
Thus  we  believe  the  assumption  of  normality  is  justified 
in  the  1968  studies  for  the  estimation  of  travel  times. 

This  should  also  be  true  for  the  estimation  of  station 
corrections  for  those  stations  that  had  large  numbers  of 
observations.  A  third  approach  is  to  employ  tests  that  are 
robust;  that  is,  perform  adequately  when  the  assumption  of 
normality  does  not  hold.  There  is  a  vast  literature  on 
the  subject  of  robustness  (Govindara julu  and  Leslie,  1970) . 
Much  of  this  work  is  on  specific  problems  and  employs  Monte 
Carlo  studies.  The  results  vary.  In  some  cases  the  usual 
test  procedures,  t,  F,  analysis  of  variance,  etc.,  work 
well  and  in  some  they  do  not.  Now  not  only  sample  size  is 
important,  but  also  population  size  since  randomization 
produces  correlation  and  the  standard  assumption  is  of 
independent  normal  observations.  Since  there  exist  broad 
classes  of  data  for  which  the  usual  tests  are  robust  (with 
suitable  sample  and  population  size)  we  again  believe  that 
the  usual  test  procedures  are  adequate  in  the  1968  studies 
in  those  cases  where  the  sample  size  was  large. 


I 
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To  compliment  these  considerations  histograms  of  the 
travel-time  residuals  were  plotted  for  each  distance  cell 
and  for  ten  selected  stations.  None  of  the  plots  exhibited 
extreme  non-normality.  We  made  no  goodness-of-f it  tests 
since  these  must  necessarily  be  chi-square  tests  (popula¬ 
tion  parameters  were  not  known). 

With  respect  to  the  Veith-Clawson  data,  we  can  again 
expect  the  usual  tests  to  be  robust  since  the  fewest  num¬ 
ber  of  observations,  50,  occurred  in  the  95°-100°  cell. 

Thus  the  t-test  indicates  that  at  a  significance  level  of 
.01  that  the  mean  of  cell  35°-40°  is  not  zero.  However, 
as  stated  earlier,  any  travel-time  function  must  possess 
certain  differential  properties,  and  the  t-test  does  not 
indicate  how  this  can  be  done.  Further,  for  these  data, 
the  t-test  is  not  exactly  correct  since  the  data  are 
corrected  residuals  and  hence  correlated  at  each  distance, 1. 
Since  optimum  estimators  can  be  obtained  without  the 
assumption  of  normality,  we  shall  do  this  and,  in  any 
event,  we  will  have  valid  estimates. 

We  suppose  that  the  Veith-Clawson  data  constitute  a 
random  sample  from  a  suitable  population.  Due  to  the  path 
averaging  this  assumption  should  hold  approximately.  The 


$ 
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fact  that  no  Russian  or  Chinese  stations  are  represented 
in  the  data,  and  all  events  are  in  the  northern  hemisphere 
implies  that  the  population  is  even  more  restricted  than 
that  of  the  .1968  studies.  These  data  represent  the  most 
extensive  collection  of  teleseismic  shot  data’ ana lyzed . 

In  particular  it  includes  all  shots  employed  by  Lilwall 
and  Douglas  (1970),  Cleary  and  Muirhead  (1969),  and 


Muirhead  and  Cleary  (1969).  These  works  are  especially  impor 
tant  since  the  results  indicate  that  the  1968  times  are  in 


error.  Ignoring  a  baseline  shift  the  corrections  with 


respect  to  these  works  are  (Douglas,  1970): 

TlDcT 1363  “  0.004* 

and 


where  LD  indicates  Lilwall-Douglas  and  CM  indicates  Cleary- 
Muirhead  and  T  denotes  travel  time. 

We  shall  follow  these  authors  and  suppose  that  any 
error  in  the  1968  times  is  given  by  a  linear  function,  say, 
a  +  In  order  to  make  comparisons  with  the  Cleary- 

Muirhead  results,  we  shall  consider  the  distance  range 
20  -100  .  Also,  we  must  allow  for  a  mean  source  correction. 


Our  model  is  then 
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=  Oj  +  T(^j)  4  a.  +  t ^  ^  4 

e//a) 

where  ^  ^  . 

T  (A^j)  =  t(  A^j)  +  a 4  bA^ 

T  (  Am)  =  ibc  I  3 68  P  i ravel  4iWj 

and  _  q  source,  ■fimc  c<?r r£c+f’oH +er^. 


Let 


Then  equation  (44)  becomes 


Uj  of  4  T  (&Q)>  4  a  4  bA^  + 


(45) 


Denote  an  observed  residual  by  Yij  (j^)  ;  that  is, 

#j ft)  =  a  +  b  Af  +  jzj  4  e*j  (lj  ,  (46) 

Equation  (46)  gives  a  suitable  model  for  estimating  a 
"smooth"  change  in  the  1968  times  by  employing  the  Veith- 


Clawson  data. 

With  randomization  our  previous  work  indicates  that 
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least  squares  analysis  is  appropriate.  Since 

+  e*j0i 

where 

g.j(j  =  2.y;jW/nj  i 

^,j)=  ;  &tu/.)=& t 

nri  ^  V 


n j  =  nuinber  of  observations  on  the  jtn  event, 


we  obtain  the  result 


-  y. j(9  a  t  < ArA-^  f  feJjT  e*jY-)). 


From  equation  (47)  we  have  directly  that 


b  =  [  •tiOl  ~"ZVvjj) 

h  (A*u'p- A'^ 

_  &  -yfyA/uj) 

j))Z 

We  are  now  confronted  with  a  dilemma.  It  is  not  possible 
to  evaluate  the  denominator  of  equat.ion  (48)  from  the  Veit’n- 
Clawson  data.  Since 

-  j)  0  ) 

(as  can  be  readily  shown),  we  have  that 


k  =  the  number  of  events. 

We  obtain  the  result  that 

b  =  0.0013 

and,  therefore, 

b  £  O.OOI3, 

a  result  in  direct  opposition  to  the  results  of  Lilwall- 
Douglas  and  Cleary-Muirhead . 

Since  inspection  of  Figure  1  indicates  that  the 
variance  of  (jJ)  decreases  with  increasing  # ,  employing 

weighted  least  squares  would  seem  advisable.  We  employ 
weights  proportional  to  the  cell  variances  to  obtain 


(49) 
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where 


and 


'/=  Var(t;,j  ft)). 


(i,j)  indicates  a  sum  over  all  i, j  at  distance  l.  It 
is  now  impossible  to  evaluate  both  the  numerator 
and  denominator  of  equation  (49)  from  the  Veith-Clawson 
data.  In  order  to  proceed  further,  we  suppose  that 

a  _  £  A*/c/ 

m_  1___ — _ — L 

J  £  Va-2 

/<5X 

(50, 

2 

Since  the  are  unknown,  we  shall  estimate  them  with 
the  data  of  Table  2.  With  these  assumptions,  we  obtain 


and 


A.=  £4.4 

t  - -0-0032 . 


However,  most  of  the  events  had  stations  over  a  wide  dis¬ 
tance  range  so  that  the  assumption  may  not 

be  seriously  in  error.  Unfortunately,  it  is  impossible 
to  ascertain  the  effect  of  this  assumption  in  equation  (49) 
Since  the  cell  residuals  are  correlated,  we  cannot 


:mmp 


even  ©S£  ipafee  r wa  H  a nee§  © f  p  ©r  b.  Thus  we  cann&t  P?ahe 
a«y  tests  even  ii  ih©  assuKptieo  ©ff  nsFsrality  is  pade.  what 
is  needed  is  the  veith*€law@©n  daft©  without  the  source  term 
removed .  uofertviately*  ikNa  data  were  unavailable.  la  (he 
afe>;i  section  we  shall  mohe  an  spproxlmatl4A  that  allows  §ta* 
list teal  tests  t©  he  mad®*  r©r  the  present,,  quel  ((alive 
conclusions  m®  be  made.  There  is  almost  certainly  a  decrease 
in  variance  with  increasing  distance.  Thus  p@r@  preei.se 
est  imators  are  obtained  by  employing  weighted  least  impure®. 

$  difficulty  in  applying  weighted  least  squares  is  that  (he 
weights  are  net  fcnewoo  ft  is  mt  clear  h§w  using  estimated 
weight®  affect®  the  respite*  Allowing  for  the  approximations 
employed  the  eel  Smite  ©f  the  slope  Is  ©.@@13  wit  hast  weight® 
and  ~@.©®3£  with  weights.  Iff  it  were  mt  fer  correlation  the 
estimated  weights  would  he  acceptable  since  the  fewest  ©b* 
serva liens  In  a  cell  ^s  1®,  in  view  ©f  these  considerations 
we  believe  that  (He  slope  is  snail  and  negative,  further, 
the  estimate  *©.@©31  is  probably  significant  or  Just  sign!® 
(leant  in  view  ©f  the  results  of  the  cell  t® tests.  $©ie 
also  that  if  *0.0031  Is  a  significant  difference  then  the 
LilwalL-oeuglas  value  off  ®0.@@6  is  significantly  different 
frow  the  value  indicated  by  the  Veiih-Clawsen  data  and  the 
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o 


o 


valu@  ©f  -@«§§3  a§  given  by  eksfy*Myifhe§^  i§  unquestionably 
significant*  that  is*  @n  a  statistical  basis  It  is  ipp@s§ible 
t©  accept  ttee  value  -@»©23  as  compared  t©  ser@9  m  eenslui© 
tbit  there  is  p@§sitely  a  slight  negatively  tilted  correction 


t©  tb®  1968  P  tip?©  will)  respect  t©  til® 
tb®  v@ltte*€lsaMs©n  Silvio  oewevtr*  w@  d@  mi  believe  t 


eorreetlsn  is  even  is  large  is  magnitude  as  Souglas  {1979)  sufoestso 
we  shill  return  to  ibis  point  is  ibe  nest  §eeii@no 


®io®t  ttee  C leary-sto i rte§ad  fioofl  results  are  based  @n 


tonisteet  residuals  ii  is  now  apparent  that  if  ttee  l%t 


times  were  slitfsl  t©  fit  the  iongoteot  data*  then  they 


e«uld  not  p@ssi.tely  fit  tbs  Veiite*cl9M«on  data.  m  fspba- 
sis©  the  fact  tbit  bow  well  a  given  world-wide  average 
travel  time  function  legates  a  specific  event  Is  not  a 
noi table  criterion  ©f  the  validity  ©f  the  travel- tin?© 
function.  in  order  to  locate  accurately  a«d  precisely  the 
effect  ©f  the  Kfi<S  }*  terw  Mist  fee  accounted  for. 

Witte  ttee  assumed  models  and  randedfca&ien  scheme  we 
have  «w©  possible  extreme  eases*  Cl)  a  simplified  model 
$®es*  in  fact*  agree  with  real  data  and  measurement  error© 
are  the  only  errors  and  fll)  there  Is  ns  measurement  error 
and  all  error©  are  due  i©  inadequacies  of  a  specific  model. 
The  actual  situation  is  almost  certainly  between  these  two 


o 


o 


*§4* 


ex tr ewes.  If  we  asnooe  rase  (i)  and  a  set  of  stat lens  fined 

^  w  t©  ©  least  squares  location 
procedure  in  ti m  errors  (Tueker,  et  alos  ig§§)0  typi^uy 

iw  le€@u@®  p&blem  the  data  ®npi@^d  arises  frrai  sta« 
ii@^§  Mast,  actually  efeserve  iht  event,  if  this  s$  €#8e> 

tfm*  it  ean*®t  fee  elated  that  the  stating  w#fe  fined  ard 
eh§e^  feef@re  observation.  However*  ear  results  indicate  that 
with  au  adequate  total  sa^l#  else  the  ordinary  est  footers  tray 
still  fee  e«?pl@yfdo  ®n  this  feasis,  the  le-eat  to  procedure  {Tusker* 
ft  al.«  19(9)  rapioyed  in  the  19i«  studies  ea®  fee  Justified. 

The  Velth-clouton  data  indicate  that  for  the  larger 
sS^ts  and  earthquakes  that  ease  fli)  mf  fee  *we  appropriate. 

If  we  assuee  ease  fl£|  we  are  led  t©  a  least  squares  leca- 
tieu  proeedure  in  dista»«ae  errors .  iff  w@  consider  the 
distauee  range  25-  *  95®  In  figure  1,  we  observe  that 
the  smoothed,  standard  deviation  is  very  nearly  proparfcion- 
al  t©  the  travel  tirre  derivative  curve*  Thus,  If  we 
divide  by  the  travel  tie®  derivative,  then  we  can 
essentially  obtain  a  constant  variance  for  the  error 
term*  that  is,  we  perform  a  variance  stabilising  transforma- 
U°n*  Trans  formations  of  this  type  are  frequently  employed 
anti  can  lead  to  Improved  results  (Kcmpthorne,  1952). 


The  fact  than  the  standard  deviation  is  pr©p®rti©nal  t© 

Hi®  travel  tin®  derivative  s§  net  nnF©a§©n#hlee  •test  iffa)* 
fmm  should  m  Mil,  If  this  is  Hi®  ease*  then  H  ffc)* 
ea©  I)®  reflected  is  a  ©mall  ditianee  change,  i&nfshot  is 
fsossstily  ns  estrone  ©a§#  and  ihe  distance  ficer  is  ehnac 
2S  Iip.  &  spall  dietenee  ehanf®  eels  like  a  per  leftist  ien 

o 

#©$  produces  a  line  ehanf#  at  a  fives  distance  proportional 
I©  the  derivative,  if  we  take  a  c®l lection  @f  these  par* 
tarhed  line  leris©  Uses  itself  standard  deviation  win  ale© 


further,  consideration  ©f  the  standard  location 


equations  (fucker,  ei  at*  19681  shews  that  dividing  hy  the 
lis»  derivative  essentially  produces  ©  location  procedure 


I©  distance.  Per  simplicity  we  ^ssuee  that  depth  is  Known 
amt  only  desire  an  epicenter  determination. (& 
an  observed  distance  and  an  initial  distance.  If 


(Tucker,  et  al< 

5?s 

$0&  8  & 

Ho 


the  true  origin  time 
=  east  longitude  of  the  event 


$  *  seis«@l©fleal  ©©-latiiutl©  ©f  the  event* 
then  ©n@  readily  Plains  the  first  @rder  result  that 

*S  &  *Jf  Sfr+(-si’'ty?+((0*Z)**+*j 

where 

6  s  the  ©rre»r  letw  in  distance 
and 

s  «■  the  aeiirosth  angle  tf-om  the  neridlan  thrnegh 
the  epicenter*  Measured  north  through  east* 
t@  the  are  free  epicenter  t©  the  station, 
equation  (11)  is  anal@'t@u§  4©  the  tine  equation  of 
Tocher*  ©t  al  (1969}  • 

Since  T(A)  is  Monotone  increasing  m3  possesses  a 
derivative  for  ©very  distance »  w©  have  that 


and 

$$At-/«.-T&g. 

Thun,  if  we  locate  by  iteration,  depth  in  Known  and  a 


common  net  of  stations  is  employed,  then  location  by 
distance  is  essentially  equivalent  to  location  by  time. 


»5?» 


fli§f§fer§,  if  distance  procedure  had  been  employed  in 
the  19&§  studies,  essentially  the  §ap§  eetiwaieg  would 
Have  been  obtained,  Slight  difference?  in  l@eati@n  would 
probably  Have  eceurred  sine#  truncal ien  w@s  e*pi@yed  In 
the  196®  studies. 

frequently  depth  is  mi  %mm  and  be 
esiifrated  and  In  reetin®  I©?#! leu  procedure®  sruncatten 
putt  always  be  eanpleytd  m  that  these  tv®  tecHniquea, 
in  general*  yield  different  leeatime  f©r  identical  initial 
date  sets.  Inspection  of  the  VeiiH-€law?®n  data  indicate? 
that  for  lit©  §H@rier  distance  it  Is  passible  le  truncate 
©  station  due  t©  podel  inadequacies  a§  well  Tg  tine  con> 
lamination,  VeliH  (personal  eoaponlcafcion,,  1.970)  report® 
that  frequently  In  routine  locations  the  closer  stations 
are.  in  fact,  truncated.  IMe  of  the  distance  location  pro* 
eedure  has  improved  many  of  these  locations  (Vcith,  per¬ 
sonal  cc*©unication,  1970).  Therefore,  it  appears  that 
the  distance  procedure  may  be  superior  to  the  tine  pro¬ 
cedure  in  large  scale  routine  location  techniques. 

Regardless  of  the  final  estimation  scheme  employed, 
the  presence  of  non-nonnol  errors  (gross  timing  errors, 
etc.)  must  somehow  bo  taken  into  account.  Thus,  wo  trun¬ 
cate  any  data  set  using  essentially  the  sample  median  and 
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$evie<ti@n  ePeet  the  j»e$ien.  We  refer  IHe  reefer 
f@  Herrin,  et  el  (19M),  page  1282  f^-c  fcHe  detail®,  tl»§ 
ipp^rtant  p@int  t@  nete  is  Uiat  any  traneatien  precede?® 
e#n  treneate  dee  t@  ©@niapi®aii@®  and  als®  dee  i®  pedel 
misfit*  In  view  @f  sHe  data  §ti  eppleyed  In  it*®  1888 
studies  we  i^etiev®  tHet  ifee  truncal i@n  da®  t@  weilel  misfit 
esppfised  a  spell  p@fti@n  ©f  tne  tfuncsted  ©teservatiens* 
§i9if@  ii§©  ^i.sjsMe  l@eati®»  pf@ged»f@  is  essentially  a 
variance  stsfoilisinf  teetsni^©®*  if  stomild  reduce  even 
fwtf^er  tls®  presort  ton  ®f  irswcailon  da®  t@  model  misfit* 


KSTIItftTIOK  ©F  U©PI»D=MIOK  TRAVIS  TJTHBS 


me  ©F  FIRST  OROBR  station  iHTRimeTiets 

F^f€TI@?i 

Assuming  (list  (lie  only  §rr§r§  are  mm:uemnt  and 

(hose  (©  Fandsnisatien  ever  a  source  effeel  and  tliac 

a  eel  of  stations  was  fined  feefere  eHservatien  {TtoeKer, 

et  al#  19691,  m  have  estimated  eerreetiens  te  (lie  Jeffreys* 

fallen  travel  tisres  (Herrin,  et  al,  1.9691  and  a  first 

erder  station  Interaction  function  fHerrle  and  Taggart, 

1969a}.  sine#  not  all  slat  lens  see  all  events,  m  Hell eve 

(Hat  (He  slat  lens  sHoald  8*e  considered  as  Having  fceen  selected 
at  randan 

him  a  fined  total  set.  0  ■••«««••  ®f  randenlsa* 

lien  seHtme  Is  ass^@d,  @wr  present  resalts  indicate  tHat 

tHe  estimate®  are  essentially  unbiased  and  optimal  Iff  m 
@§gi®@  (Hat  tlie  events  Have  Hnewn  locations.  As  previously 
observed,  Iff  we  assuage  normality,  (Hen  tHe  usual  tests 
can  He  reads. 

Since  most  ef  (He  data  in  (He  19&9  studies  were 
frens  ear tHtjualf.es ,  (He  osodel  employed  must  Include  location 

o 

terns.  How  the  first  order  agigsutHal  source  teres  are 
(equation  (21)  and  (22))* 
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snd  th#  first  order  epieenter  l©eat-l©h  terms  are  (Toeker, 
et  Hi,  1060)! 

<^9  tint  t@  S  first  ©F$§f  tin  i@gaii@n  @«$  §®yfge  effects 

e@hf@md®t?*  Is  fsefc*  if  si  I  stations  ©gcwt  st  tin 
tiisianee  from  an  events  then  the  terme  ere  pr@p@rti®nal 
the  watrl*  ©f  n@twal  e$ti aliens  feqvatiens  ©f  coed  it  lent 
i&  sififvisr.  in  this  ease,  the  effeet©  are  eespleteiy 
c@nf©3r<n$g$  @n$  ©s  estimate  @f  eaeh  separately  ean  In  stii#&0 
I©  the  19Sfe  studies,  si  I  ©vest©  employed  mere  l@©at©$ 
shots  whose  ig^ati^g  mere  knew©!  §@  tint  there  my 
h»  hia§  te  mislocetien  in  the  eet totes  ehtsineci. 

&**  previous  analysis  ©f  the  Velth-Clawson  data  indicates 
thit  for  the  travel® time  estimates  this  is  probably  ns>t 
th®  esse.  We  continue  e^r  invest  if  alien  ©f  passible 
hies  in  tin  t revel* time  estimates •  w@  shell  ignore  the 
source  mm  remove!  in  the  Veith-€laweon  data.  The  pBdel 
becomes 

3jksMp^  +  e*ic 


or 
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Where  vr-  •-3'P§y‘»<e  t|tf$ Jt  rite  are  independent  and  for  test 


purposes  nerval  ly  distributed,  Also,  equation  (5?)  1§  in 
th®  f@FP  @f  standard  regression  analysis  and  <1?©  (weighted!} 

least.  squares  estivates  are 

(jj?(  MMSS 

and  - 

Hies*  i§  I®  equivalent  t#  b*  a  fact  that  is  obvious  if 


we  ©onsider  the  least  squares 


$s»  It  Is  nsw 


is  estivate 

<* 


i 

if  *  ib®  variance  sf  p  ,  we 


ff|  */.5a36X/0'6 

and  *"  2* 

6*  sj,+/4xta' 

Mere  we  sssuve  the  waists  are  hnoMn  and  are  equal  t@  the 
estivated  eel!  variances.  Hie  observed  s  is 


%  ©  -0. 


If  we  choose  «C  ©  .01 


critical  s  is  2. 


a  iw@> 


the  critical  s  is  1. 


Hies  with  an  *  of  .01  we  do  not  reject  the  hypothesis  that 


ane<  @f  .©I  wa  do  reject.  It 


P  is  mm  and  with  an«C  @f  .©S  we  do  reject.  It  appear® 

that  there  awty  he  a  slight  correct  ten  needed  to  the  106© 

tines  with  respect  to  the  population  represented  by  the  Veltlr 
Clawson  data  or 

T*  *  fi«  -  °-eeja  * 


**  a  curate  estimate  ©f  a  teles©  iiu®Je  shot  world-wide 
average 

4travel-tim@  function.  ll^ever,  the  results  dre  not  that 
conclusive  end  we  Relieve  that*  tossed  on  ttoe  Veith-Clawson 
dsts,  ttoere  Is  no  compelling  reason  to  slier  ttoe  1960 
iisres.  Note  sis©  that  If  we  test  the  hypothesis  that 0 
is  ser©  versus  the  alternative  that^  5  -0.006  (the  L-D 
value!  ve  obtain  analogous  results.  ?he  observed  e  is 
again  -2.24  tout  now  the  test  is  one-sided  so  that  the 
critical  »* s  are  -2.326  (0C  =  .011  and  -1.645  (o<  s  .09). 

Por  comparison,  we  shall  test 

il0*  p  *  -0.005  (L-D) 
versus 

M*«  £  >-0.006 
and 

II  *  |9  s  “0.023  (C-M) 
versus 

llAi  I?  >-0.023. 

Por  the  first  test  the  observed  s  is  1.98.  the  critical 
8*®  8re  2*236  1*=  •01>  and  1.W5  («s  ,o§).  (7|te  test 
is  in  the  positive  direction.)  We  again  reach  the  same 
conclusion  with  respect  to  the  value  suggested  toy  Lilwall 
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and  Douglas  (Douglas,  1970).  However,  f or  the  second  tost 
the  observed  ?.  is  14.00  which  is  significant  with  any 
reasonable  o< level.  The  value  of  -0.023  suggested  by  Muir- 
head  and  Cleary  (1969)  is  just  not  consistent  with  the 
Vcith-Clawson  data. 

To  further  test  for  possible  bias  due  to  dislocation 
in  the  travel  time  estimates  we  made  n"Monte  Carlo"  study 
of  island  arc  events  in  the  Pacific  and  Atlantic.  The  actual 
events  simulated  and  the  procedure  is  based  on  a  model  of 
Longshot  which  will  be  discussed  in  detail  in  the  section 
on  the  source  interaction  function.  The  event  regions  and 
events  for  this  study  are  given  in  Table  3.  We  employed 
20  events  and  52  stations.  The  stations  are  given  in  Table  4 
All  stations  employed  are  from  the  1968  studies  (Herrin  and 
Taggart,  1968a) .  We  attempted  to  make  the  station  selection 
both  world-wide  and  representative  of  those  stations  for 
which  we  have  estimated  station  corrections.  The  event 
regions  are  basically  island  arc  regions  and  all  but  two 
ore  in  the  Pacific.  We  attempted  to  make  these  regions 
representative  of  the  data  set  employed  in  the  1968  studies. 
Most  of  the  event  regions  are  in  the  western  pacific  which 
is  consistent  with  the  1968  data  set  (Herrin,  et  al.,  1968). 
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The  Puerto  Rienn  event  in  represnetative  of  only  one  event 
near  the  Lessor  Antilles  and  as  such  might  not  hove  boon 
included.  To  the  extent  that  these  events  and  stations  are 
representative  of  the  1968  data  the  results  will  indicate 
the  possible  effects  of  mislocation  in  travel  time  bias.  Wc 
tacitly  assumed  that  event  locations  outside  of  island  arc 
regions  in  the  1968  studies  were  unbiased  and  hence  no  data 
for  these  regions  is  included  in  the  study.  To  give  an  ex¬ 
treme  case  we  assumed  that  a  source  function  equivalent  to 
Longshot  existed  in  each  event  region.  This  source  function 
approximates  that  of  the  Longshot  model  study  and  is  given  by 


[-0.5  +  |  .32  sin  +275.6°+  6>k)]  (Herrin  and  Taggart, 
1968b)  where  0^  gives  an  adjustment  to  for  the  trend  of 

the  island  arc.  (For  the  Longshot  event  0#  )  With 

this  source  function  and  the  1968  times  we  generated  observed 
arrival  times  at  each  of  the  stations.  No  additional  random 
errors  were  introduced  so  that  the  study  is  Monte  Carlo  in 
nature  only  to  the  extent  that  the  event  regions  and  stations 
constitute  a  random  sample  from  the  1968  data  set. 

The  generated  arrival  times  were  then  employed  in  .locat¬ 
ing  the  events.  Time  residuals  were  obtained  from  these 
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estimated  locations.  The  residuals  were  then  grouped  on  a 
5°  cell  basis  to  compute  the  mean  and  otondard  deviation  and 
the  results  arc  presented  in  Table  5.  Tho  means  given  in 
Table  5  can  be  thought  of  as  estimated  corrections  to  the 
i960  times.  We  thus  tafcc  as  true  the  1968  times  and  use 
data  that  has  as  error  only  model  misfit  cr  source  interaction. 
We  then  simulated  the  first  iteration  in  a  procedure  to  estimate 
the  true  travel  times.  This  should  represent  an  extreme  case 
since  only  one  step  in  the  iteration  was  performed.  In  find¬ 
ing  the  final  estimated  curve  we  would  apploy  the  corrections 
given  in  Table  5  to  .obtain  a  new  travel-time  curve,  then  re- 
estimate  the  locations,  etc.  The  data  of  Table  5  represent 
that  portion  of  estimated  corrections  that  would  be  due  to 
mis location  and  as  such  are  over-corrections.  Any  estimates 
that  have  this  phenomenon,  included  should  have  the  corrections 
of  Table  5  subtracted  in  order  to  eliminate  mis  location  bias. 

The  results  of  Table  5  indicate  that  the  estimated  curve 
would  be  slow  with  respect  to  the  true  curve  from  20°  to  60°, 
unchanged  in  the  range  60° -75°  and  fast  in  the  range  75° -100  . 
The  estimated  slope  using  weighted  least  squares  with  the 
observed  variances  as  weights  yields  a  slope  estimate  of 
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-0.0030.  However,  this  is  a  positive  slope  with  respect  to 
the  Vcith-Clawson,  Lilwall-Douglas  or  ^leary-Muirhead  data 
and  indicates  that  any  rotation  produced  in  the  1968  times 
from  this  effect  is  in  a  direction  opposite  to  that  suggested 
by  Douglas  (1970)  and  Muirhead  and  Cleary  (1963).  While  only 
one  iteration  was  made,  it  is  unlikely  that  further  iterations 
would  change  the  slope  from  +0.0030  to  -0.0032  or  more  negative 
values.  If  the  Longshot  model  represents  an  extreme  case  of 
source  function  and  the  selected  events  are  representative 
of  the  1968  study  then  we  conclude  that  any  bias  due  to  mis- 
location  should  be  small.  The  small  negative  slope  indicated 
by  the  Veith-Clawson  data  could  be  due  to  effects  other  than 
mis location. 

As  a  final  check  we  can  compare  the  1968  times  to  the 

Lilwall  and  Douglas  (1970)  times.  Employing  the  smoothed 

times  and  95%  confidence  intervals  as  given  by  Lilwall  and 

Douglas  (1970)  for  their  estimated  travel-time  curve  and 

the  1968  times  and  95%  confidence  intervals  (Herrin,  et  al., 

1968),  one  can  perform  a  significance  test.  Ignoring  a 

baseline  shift  the  confidence  bands  of  the  two  curves 

(1970)  curve 

(Lilwall  and  DouglasAconsidered  at  only  their  smoothed 
points)  fail  to  touch  or  overlap  at  only  the  distance  43°. 


f. 
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The  separation  is  approximately  0.1  sec.  Since  adjacent 
confidence  intervals  are  not  independent,  it  is  impossible 
to  determine  the  significance  level  of  the  test.  Further, 
the  data  sets  overlap  somewhat,  and  many  criteria  present 
themselves.  We  could  reject  the  hypothesis  that  the  two 
curves  represented  by  the  two  data  sets  are  equivalent  if 
one  or  more  confidence  intervals  do  not  touch  or  overlap  or 
if  two  or  more  do  not  touch  or  overlap,  etc.  It  is  not 
clear  which  criterion  is  appropriate.  We  submit  that  the 
data  do  not  suggest  that  the  two  curves  differ  seriously. 

Lilwall  and  Douglas  (1970)  and  Douglas  (1970)  claim 
that  the  rotation  between  the  L“D  times  and  the  1968  times 
is  possibly  due  to  event  mislocation.  In  support  of  this 
claim,  they  simulated  one  iteration  of  the  1968  studies  with 
their  data.  However,  it  is  not  correct  to  compare  on  just 
one  iteration.  One  iteration  gives  estimates  of  the  travel 
times  that  are  unadjusted  for  location  and  other  model  effects 
Any  unbiased  estimates  must  be  adjusted  for  all  effects. 

This  is  analogous  to  attempting  to  invert  a  matrix  and  not 
accounting  for  the  off  diagonal  terms.  The  total  estimation 
scheme  employed  in  the  1968  studies  (Tucker,  et  al.,  1968) 
was  essentially  a  Siedel  iterative  procedure  (Householder, 
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1953).  The  siede.l  procedure  will  converge;  that  is,  give 
the  same  estimates  as  matrix  inversion  if  certain  conditions 
are  met.  We  believe  these  conditions  were  met  in  the  1968 
studies  and  probably  for  the  Lilwall-Douglas  (1970)  data. 

If  this  is  the  case,  then  the  joint  epicenter  method  of  Lil- 
wall  and  Douglas  and  the  Siedel  method  of  the  1968  studies 
are  mathematically  equivalent.  In  order  to  compare  correctly 
the  two  methods,  the  Lilwall-Douglas  data  should  have  been 
iterated  enough  times  to  be  equivalent  to  the  twelve  itera¬ 
tions  performed  in  the  1968  studies. 

We  conclude  that  there  is  no  appreciable  bias  in  the 
1968  times  and,  hence,  at  present  no  compelling  reason  to 
alter  these  times.  As  more  shot  data  or  data  from  large, 
well  located  earthquakes  becomes  available,  it  should  be 
possible  to  estimate  any  change  in  the  1968  times.  Since 
it  appears  that  a  slight  rotation  may  be  needed  to  make  the 
1968  times  more  nearly  fit  a  teleseismic  shot  average, 
we  offer  a  possible  explanation.  The  indicated  change  is 
probably  due  to  many  sources.  We  first  believe  that  there 
may  be  rea!l  differences  between  a  travel-time  curve  for  the 
population  represented  by  the  Veith-Clawson  shot  data  and 
the  population  represented  by  the  1968  data.  Inspection  of 
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Table  1  indicates  that  only  Longshot  or  approximately  5%  of 

the  data  come  from  active  seismic  regions.  Thus  the  Veith- 

Clawson  population  is  primarily  from  aseismic  areas  while  that  of  th 

earthquake  zones. 

1968  studies  is  primarily  from  /  The  important  result,  we 
believe,  is  that  the  Ve ith-clawson  data  indicate  that  the 
1968  times,  based  essentially  on  earthquake  data,  are  also 
applicable  to  shots  in  aseismic  regions.  Also  there  is 
possibly  a  slight  rotation  due  to  mislocation.  However,  part 
of  any  rotational  change  may  be  due  to  the  arbitrary  model 
employed  for  the  distance  range  0°-20°  in  the  1968  studies. 
Inspection  of  the  distance  ranges  10° -15°  and  l5°-20°  of  the 
Ve ith-clawson  data  indicates  that  the  1968  times  are  highly 
significant  and  fast.  Since  a  change  in  the  model  could  pro¬ 
duce  a  slight  rotation  (Herrin,  et  al. ,  1968)  and  the  effects 
of  the  arbitrary  model  may  go  out  to  30°,  a  different  initial 
model  can  account  for  part  of  any  rotational  change.  Finally, 
regional  variations  in  P  times  may  occur  at  distances  as 
great  as  60°.  If  this  is  the  case,  then  estimates  in  the 
range  25°-60°  could  vary  widely  depending  on  the  actual  data 
set . 

While  there  is  apparently  only  slight  or  no  bias  in  the 
1968  times,  the  results  on  station  terms  are  not  as  conclusive. 
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Witli  the  exception  of  the  Lilwall-Douglas  (1970)  work,  the 
station  corrections  show  the  same  overall  features  as  other 
published  corrections.  The  Lilwall-Douglas  corrections  for 
European  stations  appear  to  have  about  a  0.5  sec.  negative 
displacement  with  respect  to  those  of  the  1968  studies.  This 
shift  may  or  may  not  be  statistically  significant.  Further, 
the  randomization  analysis  indicates  that  a  minimum  of  50  to 
100  observations  per  station  are  required  in  order  to  have 
the  order  of  the  bias  be  0.02  sec  or  less .  We  believe  that 
those  station  terms  of  the  1968  studies  based  on  50  or 

more  observations  are  adequate.  We  reason  as  follows:  To 
the  extent  that  there  is  no  travel  time  bias,  the  location 
terms  were  completly  confounded  with  any  first  order  azimuthal 
source  terms  and  hence  were  actually  estimated.  We  cannot 
separate  a  mislocation  and  an  azimuthal  source  effect.  Thus, 
any  estimable  model  will  not  have  separate  location  and  azi¬ 
muthal  source  terms,  one  or  the  other  or  some  linear  com¬ 
bination  of  the  two  will  suffice.  Since  a  mislocation  accounts 
for  nearly  all  of  a  first  order  azimuthal  source  effect,  there 
are  no  bias  terms  in  the  model;  that  is,  to  a  first  order  a 
possible  adequate  model  has  (i)  location  terms,  (ii)  station 
effect  terms  and  (iii)  possibly  travel  time  terms. 


PROBLEM  OF  THE  SOURCE  INTERACTION  FUNCTION 


Since  the  location  terms  are  confounded  with  the  first  order 
azimuthal  terms  of  a  source  interaction  function,  in  general,  we 
mis located  in  the  presence  of  strong  source  terms.  While  this 
was  desirable  for  the  1968  studies  it  is  undesirable  for  a  pure 
location  problem.  In  order  to  investigate  the  possibility  of 
estimating  a  source  function  and  a  location  simultaneously  we 
performed  a  Monte  Carlo  study  using  the  estimated  Longshot 
source  function  (Herrin  and  Taggart,  1968b)  and  a  network  of 
stations  from  those  employed  in  the  1968  studies.  The  network 
was  balanced  to  the  extent  possible  and  had  stations  well  dis¬ 
tributed  in  azimuth  and  distance  from  the  source.  We  assumed 
an  event  located  on  Amchitka .  No  station  terms  were  included 
or  estimated.  The  46  stations  employed  are  given  in  Table  6. 

The  study  was  performed  as  follows:  With  the  fixed  network 
the  Longshot  model  was  used  to  derive  time  terms  for  each  station. 
These  were  added  to  the  true  travel  times  (as  given  by  the  1968 
P  times)  from  event  to  station.  This  constituted  a  fixed  set 
of  initial  arrival  times  for  each  station.  To  these  times  were 
added  independent  normal  deviations  with  constant  variance  O* , 
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for  all  stations.  The  values  of  (5  =0,  0.08,  0.8  and  1.0  sec 
were  employed.  The  initial  arrival  times  plus  the  normal  error 
terms  were  employed  as  a  set  of  observed  arrival  times  at  each 
station.  These  observed  arrival  times  were  employed  in  a 
joint  location  and  source  term  estimation  procedure.  The  actual 
location  was  recorded  and  a  series  of  runs  made  for  each 
variance  level,  in  the  case  of  $  =0  only  one  run  was  made. 

The  location  was  0.065  kms  north  and  0.005  kms  west  of  the  true 
location,  that  is,  with  no  error  we  can  locate  and  estimate  a 
source  function  simultaneously.  The  matrix  of  normal  equations 
is  not  singular  in  a  neighborhood  of  the  true  location.  For 
the  other  variance  levels  repeated  runs  were  made.  The  mean 
location  change  and  standard  deviation  (both  in  kms)  were  ob¬ 
tained  for  the  north-south  and  east-west  directions  about  the 
true  location.  The  results  are  presented  in  Table  7  where 

6 N-S  “  /<£*  anc*  Gff-W  S  • 

To  check  the  validity  of  the  results . certain  tests  were 
made.  Since  multiple  minima  may  occur  in  nonlinear  least 
squares  an  extensive  search  was  made  of  a  selected  run.  For 
this  run  the  location  without  source  terms  was  16.64  km  north 
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and  9.29  lot  cost  while  the  location  with  source  terms  was 
29.81  km  south  and  5.50  ten  west.  The  location  without  source 
terms  is  apparently  typical  while  the  location  with  source 
terms  appears  to  be  too  far  south.  If  spurious  minina  are 
commonly  reached  this  run  should  serve  as  an  example.  The 
search  indicated  that  the  iterative  least  squares  procedure 
did  in  fact  reach  the  true  minimum  and  that  there  there  is 
only  one  minimum.  This  is  not  unexpected  since  the  travel 
time  nonlinearities  are  not  great.  Throughout  the  study 
different  initial  points  were  employed  on  the  same  run  in 
order  to  check  for  dependence  on  the  initial  point.  In  all 
cases  the  final  locations  were  identical.  The  more  common 
initial  point  was  the  true  location  since  this  minimized  com¬ 
puter  run  time. 

The  last  set  of  runs  given  in  Table  7  also  serve  as  a 
check.  For  this  group  no  source  function  was  estimated.  This 
results  in  locations  biased  toward  the  north.  The  "true' location, 
that  is,  the  location  without  error  is  given  by  the  run  with 
S  -  0  and  is  22.46  km  north  and  1.96  km  west.  The  mean  shifts 
of  -0.715  and  -0.179  are  not  statistically  significant  with 
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oc  a  0.05  or  cx  *0.01.  In  fact  the  N-S  shift  is  barely  sig¬ 
nificant  at  the  im  level  and  the  B-W  shift  is  not  significant 
until  the  63%  level  is  reached.  Here  we  test  each  direction 
independently.  This  is  approximately  valid  since  the  network 
of  stations  is  nearly  balanced.  The  off  diagonal  terms  in 
X'X  are  smaller  by  an  order  of  magnitude  than  the  diagonal 
terms  and  the  confidence  ellipse  is  essentially  circular, 
practically,  the  north-south  and  east-west  locations  are 
statistically  independent.  Based  on  this  consideration  the 
results  were  presented  for  the  individual  axes  in  Table  7. 

Since  these  checks  indicate  that  the  Monte  carlo  runs 
were  properly  made  we  shall  now  investigate  the  implications 
of  the  results.  First  wc  see  that  the  nonlinearities  produce 
no  apparent  effect  on  the  uriancc  of  the  estimated  locations. 
All  of  the  S  lie  essentially  between  15  and  20  km  for  those 
runs  where  source  terms  were  estimated.  Also  the  standard 
deviations  for  the  non-source  runs  were  consistent  with  real 
data  for  large  events  located  by  good  stations  with  balanced 
networks,  standard  deviations  of  5-7  km  are  commonly  observed. 
It  appears  that  introducing  source  terms  into  the  given  network 
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produced  a  fourfold  increase  in  the  standard  deviation  of  the 
estimated  location.  This  is  not  unexpected  since  even  for  a 
veil  balanced  network  the  location  terms  and  azimuthal  source 
terms  are  partially  confounded.  The  normal  equations  while  ap¬ 
parently  not  singular  have  near-singular  tendencies. 

Secondly  we  note  that  the  estimated  location  is  affected 
by  model  nonlinearities.  As  the  error  standard  deviation  is 
increased  the  bias  increases.  This  is  not  uncommon  in  non¬ 
linear  regression.  Investigation  of  the  0*0.08  case  indicates 
that  the  bias  is  not  significant  for  any  reasonable  Of  .  In 
fact  it  appears  that  for  balanced  networks  if  the  only  errors 
were  observational,  then  we  could  simultaneously  estimate  a 
location  and  source  effect. 

Finally,  with  the  balanced  network  studied  it  appears  infeas¬ 
ible  to  estimate  both  a  location  and  azimuthal  source  effect. 

The  6  *  0  case  yields  a  bias  of  22.46  km  and  a  standard 
deviation  of  3. 99  km  for  the  north-south  direction  and  a  bias 
and  standard  deviation  of  1.96  km  and  3.43  km  respectively  for 
the  east-west  direction.  The  corresponding  mean  square  errors 
are  620.36  km2  and  15.61  km2  for  the  respective  directions  or 
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an  overall  MSE  of  1335.07  km2  based  on  independence,  if  we 
take  G  «  1.0  for  a  typical  location,  then  similar  computa¬ 
tions  yield  for  the  source  estimated  case  339.59  km2  and 
294.61  km2  MSP  for  the  respective  directions  for  an  overall 
MSB  of  634.20  km2.  Thus  the  source  estimated  case  docs  not 
compare  unfavorably.  However,  let  us  now  suppose  that  there 
is  no  true  source  effect  present.  As  is  indicated  by  the  last 
case  in  Tabic  7  and  by  the  work  of  Tucker,  ot  al.,  (1968) in 
the  sourcc-not-cstimated  case  the  bias  would  bo  zero  and  the 
overall  MSB  27.56  km2.  However,  when  the  azimuthal  source  is 
estimated  the  procedure  no  longer  acts  linear.  It  is  not  clear 
what  the  bias  would  be.  A  conservative  estimate  is  zero.  Fur¬ 
ther  the  variances  should  remain  the  same  since  the  nonlinoaritic 
apparently  do  not  affect  them.  With  these  assumptions  the  MSE 
for  the  source  estimated  case  with  no  source  function  present 
is  633.89  km2.  In  this  case  it  would  be  unreasonable  to 
estimate  simultaneously  a  location  and  an  azimuthal  source  effee 
These  results  are  based  on  an  excellent  network,  prac¬ 
tical  networks  are  often  poorly  balanced.  It  is  of  interest 
to  compare  the  two  procedures  under  these  conditions.  To 
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simplify  matters  and  reduce  run  timo  only  the  zero  error  case 
was  studied.  The  networks  designated  by  I,  II ,  Til  and  IV 
and  given  in  Table  8  were  studied.  The  results  are  presented 
in  Table  9.  As  one  might  expect  the  source  estimation  pro¬ 
cedure  worked  reasonably  well  without  observation  error  (5=0). 
However,  for  network  XV  the  procedure  yielded  a  n-S  error  of 
72.93  km  as  compared  to  16.61  km  when  locating  without  source 
terms.  This  is  clearly  unacceptable.  Further,  network  IV  is 
similar  to  network  III.  In  fact  IV  was  selected  while  attempt¬ 
ing  to  obtoin  a  network  that  was  a  reflection  of  III.  There 
is  no  indication  in  the  two  networks  that  the  N-S  error  would 
go  from  0.45  km  to  72.93  km.  There  is  no  way  in  practice  to 
know  that  this  type  of  phenomenon  has  occurred. 

We  pursue  these  considerations  by  now  considering  a  real 
event,  Longshot,  which  was  large  and  well  recorded.  The  Long- 
shot  observations  were  initially  sorted  on  the  estimated  sta¬ 
tion  variances  and  the  number  of  observations  per  station  as 
given  by  Herrin,  et  a!.,  (1968).  The  criteria  were:  minimum 
number  of  observations  25  and  maximum  variance  1.00  sec2. 

Prom  this  sorting  two  reasonably  balanced  networks  were 
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s elected.  They  are  given  in  Table  10  and  are  labled  V  and  VI. 
The  fewest  observations  for  V  was  34  and  the  largest  variance 
was  1.00  sec2.  For  VI  these  quantities  are  33  and  .92  sec2 
respectively,  with  these  networks  a  number  of  combinations  of 
location  procedures  were  investigated.  For  each  network  we 
located  with  and  without  station  corrections  (Herrin,  et  al., 
1968)  and  with  and  without  estimates  of  azimuthal  source  cor¬ 
rections.  The  results  are  presented  in  terms  of  north-south 
and  east-west  errors  in  km  in  Table  11  and  are  discouraging. 

The  station  corrections  appear  to  work  reasonably  well  when  a 
source  function  is  not  estimated.  Unfortunately,  when  we  esti¬ 
mate  an  azimuthal  source  function  the  station  corrections  appear 
to  further  degrade  the  location.  However,  we  note  that  general¬ 
ly  the  corrections  are  in  the  right  directions,  even  when  all 
are  applied  or  estimated;  that  is,  the  corrections  move  the 
estimated  location  south.  Network  V  gives  a  location  19.86  km 
north  of  the  true  and  VI  gives  a  location  30.22  km  north  of 
the  true.  The  apparent  true  location  is  about  25  km  north  c£ 
the  true.  Errors  (both  measurement  and  model)  in  the  data  for 
these  two  networks  are  apparently  such  that  "precise"  solutions 
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and  Sorrells,  i969)  A  mC  re9ion  (Herrin 

-  —  -tie  materia  T;  3  d0™'90ln9  “*  —  Plate 

-epts  of  ocean-floo  ^  «- 

floor  spreading  and  global  Dlat. 

18  rePreS°nted  in  ^  by  Pigure  3.  The 

len9th  °f  «"  was  determined  by  th  3"  S  “ 

y  thS  SPat«l  distribution 
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of  earthquakes  in  the  central  Aleutian  Islands.  Using  a 
three-dimensional  ray  tracing  program  developed  by  G.  Sorrells, 
we  calculated  the  travel  time  residual  pattern  relative  to  the 
1968  travel  times  which  would  result  from  this  distribution  of 
velocities  beneath  Amchitka.  Figure  4  shows  the  residuals  as 
contours;  the  map  is  symmetrical  about  the  north— south  line 
so  that  only  the  eastern  half  is  plotted.  Note  the  complexity 
of  the  map  and  that  the  most  negative  residuals,  representing 
the  fastest  ray  paths,  are  at  distances  of  20°  to  60°  and 
azimuths  northeast  and  northwest  from  Longshot.  In  Figure  5 
the  residuals  are  plotted  as  a  function  of  distance  for  four 
azimuths.  Using  the  actual  locations  of  84  stations  which 
recorded  Longshot  we  computed  the  station  residuals  which  would 
result  directly  from  this  model.  Figure  6,  which  shows  these 
residuals  plotted  as  a  function  of  azimuth,  is  very  similar 
to  the  same  kind  of  plot  of  observed  Longshot  residuals.  A 


first  order  sine-curve  was  fit  to  the  model  residuals  of  the 
Bsin(H|V  +  €)  where  6  is  the  phase  angle. 


form 
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The  resulting  values  were 


[-0.32  +  0.6S  s.VifZj;  +  250“)] 


which  can  be  compared  to  the  fit  to  observed  Longshot  residuals 
which  was,  [“0*5+  f.3£siV2j-+  Z73.(f)]. 

These  functions  are  similar;  however,  the  observed  data  gave 
larger  values  for  both  the  A  term  and  the  amplitude  of  the  sine 
function.  The  similarity  can  be  improved  by  increasing  the 
velocity  contrast  in  the  assumed  source  model  or  by  lengthening 
the  down  dipping  plate. 

The  source  parameters  ,  Ay*  and  By  £  in  the  first 

order  model  were  calculated  for  the  source  region  represented 


by  Figure  3.  We  then  estimated  the  Longshot  epicenter  using 
equation  (20)  as  a  statistical  model,  with  the  distance  depen¬ 
dence  of  the  station  terms,  and  being 

neglected.  The  resulting  location  shown  on  Figure  7  as  a 
dot  marked  "corrected  times"  is  closer  to  the  true  location 
than  our  previous  estimate.  Had  the  length  of  the  dipping 
plate  in  our  assumed  source  model  been  300  km  instead  of  200  km, 
the  computed  location  would  have  been  very  close  to  the  true 
shot  point. 
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The  use  of  deterministic  source  models  appears  to  be 

l 

the  most  promising  course  to  follow  in  attempting  to  remove 
bias  from  epicentral  estimates  employed  in  the  1968  studies. 

Our  previous  estimates  of  azimuthally  dependent  station 
corrections  may  well  be  biased  because  of  the  spatial  corre¬ 
lation  of  source  effects  along  the  island  arcs  of  the  Pacific 
and  Atlantic.  The  next  step  in  a  re-estimation  procedure  is 
to  relocate  all  of  the  island  arc  events  we  have  used,  assum¬ 
ing  the  travel  time  tables  are  correct,  based  on  a  modified 
equation  (20)  and  assumed  source  models  for  each  arc.  (The 
source  terms  are  modeled  rather  than  estimated.)  New  residuals 
then  obtained  can  be  used  to  re-estimate  the  station  terms. 

The  process  must  involve  serveral  interations  because  the 
time  anomaly  field  for  a  given  source  model  is  very  sensitive 
to  a  shift  in  location  normal  to  the  strike  of  the  dipping  plate 


CURRENT  CAPABILITIES  IN  THE 


ESTIMATION  OF  LOCATION  PARAMETERS 


The  Monte  Carlo  study  of  Amchitka,  the  Longshot  study 
analysis  of  the  Veith-clawson  data  and  the  work  of  Tucker, 
at  al.,  (1968)  indicates  that  an  error  variance  of  1  sec2  is 
not  unreasonable  for  present  network  capabilities  and  larger 
events.  This  reflects,  primarily,  errors  due  to 
(i)  higher  order  model  inadequacies  and 

(11)  basic  station  operational  capabilities  (timing, 
reading,  etc.). 

The  overall  value  of  1  sec^  is  only  justified  on  a  suitable 
randomization  basis  since  model  inadequacies  are  not  true 
random  variables.  Further,  in  order  to  achieve  this  value 


we  must  routinely  truncate  our  data  sets  to  eliminate  con¬ 
tamination  errors,  since  the  first  error  source,  (i),  is 
approximately  an  order  of  magnitude  larger  than  the  second 
source,  (ii),  (for  the  better  stations)  any  reduction  in  ff 2 
must  come  from  more  accurate  models. 

As  a  first  step  the  distance  dependence  in  station  cor¬ 
rection  terms  should  be  accounted  for.  The  Amchitka  Monte 
Carlo  study  indicates  that  if  more  accurate  and  precise 
estimates  of  the  station  time  terms  were  available  it  might 
be  possible  to  simultaneously  locate  and  estimate  azimuthal 
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source  terns.  That  is,  if  we  can  get  the  magnitude  of  error 
source  (i)  comparable  to  that  of  source  (ii)  We  can  possibly 

eliminate  by  estimation  the  bias  due  to  an  azimuthal  source 
function.  However,  if  an  approximate,  deterministic  source 
model  is  available  (estimated  or  known  from  other  work),  then 
the  value  of  1  sec*  is  operationally  adequate  as  evidenced  by 
the  work  of  Tucker,  et  al.,  (1968),  Herrin  and  Taggart  (1968b) 
and  the  model  studies  of  our  present  work.  All  of  these 
studies  indicate  that  for  a  well  distributed  network  composed 
of  25  or  more  good  stations,  the  major  contributor  to  location 
error  is  source  bias  produced  by  an  azimuthal  source  function, 
in  the  Herrin-Taggart  (1968b)study  the  improvement  in  locations 
produced  by  use  of  station  corrections  was  of  the  order  of  a 
few  km  while  the  improvement  by  use  of  source  corrections  was 
Of  the  order  of  tens  of  km.  The  Monte  carlo  studies  of 
Tucker,  et  al.,  (i968)  gave  similar  results,  our  present 
model  studies  indicate  that  this  phenomenon  is  due  to  the 
confounding  of  the  location  terms  and  azimuthal  source  terms 
employed  in  the  least  squares  estimation  procedure.  Thus 
any  real  improvement  in  location  capabilities  must  come  from 
adequately  accounting  for  azimuthal  source  terms. 

We  have  as  direct  implications  of  the  Monte  Carlo  studv 


Of  Amclutka  and  the  Longshot  studies  that 

U)  bias  in  a  simultaneous  location  and  azimuthal 
source  term  estimation  procedure  is  of  the  order  of  a  few  km, 

(ii)  in  locating  either  with  or  without  azimuthal  source 

terms,  the  procedure  is  linear  with  respect  to  the  variance 

of  the  estimated  location  which  implies  that  the  location 

diagonal  terms  of  (X'X,-1  increase  when  estimating  azimuthal 
source  terms, 

in  either  case  a  M 

e  a  suifcable  measure  of  network 

effectiveness  is  given  by  (x'x)"1,  and 

(iv)  for  a  location  without  source  term  estimation, 

using  a  balanced  network  of  oc  _ „ 

network  of  25  or  more  good  stations,  standard 

deviations  of  5  to  7  km  can  be  obtained. 

Two  possible  solutions  to  the  source  problem  are  now 
apparent,  xn  one  we  can  use  (X'X)-l  as  a  figure  of  merit  and 
attempt  to  minimize  the  overall  MSE  by  estimating  only  a 
proportion  of  any  azimuthal  source  terms,  on  the  other  hand 
we  can  model  the  local  structure,  e.g.,  Longshot,  and  use 
this  model  to  predict  residuals.  Both  have  advantages  and 
disadvantages.  The  first  requires  no  specific  model.  »  is, 
however,  limited  to  a  few  terms  in  an  expansion  and  has  the 
drawback  of  estimating  only  a  proportion  of  any  possible 
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lower  order  azimuthal  source  terms.  This  is  partially  com¬ 
pensated  by  the  fact  that  first  order  terms  should  be 
adequate  and  a  measure  of  the  overall  MSE  is  available.  The 
proportion  of  any  source  terms  estimated  is  known  and  can  be 
used  to  predict  the  bias  and  also  the  variance  of  the  loca¬ 
tion  estimators  is  known  from  (X'X)”1  (assuming  <5  2  =  1  sec2). 
Thr  second  method  has  an  overall  variance  that  is  known  and 
adequate  in  view  of  (iv) .  Unfortunately  any  bias  removal  is 
a  direct  function  of  the  model  and  for  many  source  regions  we 
have  no  knowledge  of  proper  models.  Not  only  is  the  model 
important  but  also  the  location  of  the  event  with  respect  to 
the  model.  For  example,  with  Longshot  a  change  in  the  epicenter 
of  about  50  km  to  the  south  should  virtually  eliminate  all  of 
the  time  residuals  due  to  the  Longshot  source  function.  Also 
any  location  procedure  would  be  costly  since  a  search  tech¬ 
nique  must  be  employed.  (Analytic  representations  of  the 
models  we  have  considered  are  difficult  to  derive.) 

Either  of  these  methods  can  be  equally  effective.  Our 

Longshot  model  studies  indicate  that  a  reasonable  model 

accounted  for  about  one-half  of  the  bias.  Assuming  a  stan- 

2 

dard  deviation  of  6  km  we  obtain  an  overall  MSE  of  144  km  + 

36  km2  or  180  km2.  To  evaluate  the  first  method  let  us 
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assume  that  a  proportion  of  one-half  yields  a  standard  de¬ 
viation  of  10  to  12  km,  say  11  km.  Preliminary  work  with  the 
proportion  method  locating  the  Monte  Carlo  Amchitka  event 
indicates  that  these  assumptions  are  not  unreasonable,  with 
these  assumptions  we  obtain  an  overall  MSE  of  144  km2  +  121  km2 
or  265  km2  or  slightly  poorer  than  the  model  method.  However, 
both  MSE's  are  comparable.  Thus  we  conclude  that  with  proper 
modeling  or  careful  proportional  estimation  of  source  terms 
epicenter  location  MSE ' s  of  180  to  265  km2  are  obtainable. 

Further,  improvement  in  these  values  should  be  obtained 
by  employing  the  distance  least  squares  method.  We  have  pre¬ 
viously  discussed  the  usefulness  of  the  distance  method  in 
routine  location  procedures  where  truncation  is  employed. 

Recall  that  the  distance  method  is  essentially  the  result  of 
a  variance  stabilizing  transformation.  To  the  extent  that 
the  distance  method  produces  a  constant  variance  it  will,  in 
theory,  result  in  lower  variances  for  the  location  parameters. 
This  is  so  since  in  the  case  of  unequal  variances  a  weighted 
least  squares  is  optimum.  Thus  the  distance  method  should 
improve  least  squares  depth  estimates  by  retaining  closer 
stations  and  should  reduce  the  variance  of  all  location 


parameters . 
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Ic  is  difficult  to  assess  the  quantitative  affects  of 
the  distance  method  on  epicenter  location  errors;  however, 
the  improvement  in  standard  deviation  should  be  only  a  few 
km.  If  we  consider  the  model  case  and  suppose  that  the 
standard  deviation  reduction  is  2  km  we  obtain  144  km2  + 

16  km2  or  an  overall  MSE  of  160  km2.  Similar  considerations 
for  the  proportional  case  yield  144  km2  +  81  km2  or  an  over¬ 
all  MSE  of  225  km2.  We  then  obtain  overall  epicenter  location 
MSE ' s  of  160  to  225  km2.  Again  we  see  that  bias  reduction  is 
the  important  consideration. 

By  employing  "good"  station  corrections  a  further  re¬ 
duction  of  the  location  parameter  variances  is  possible. 

These  improvements  are  again  slight  with  respect  to  the 
possible  improvement  in  location  bias.  We  conclude  that 
by  employing  a  distance  location  procedure,  "good"  station 
corrections  and  either  a  source  model  or  the  proportional 
method  that  for  a  well  distributed  network  composed  of  25  or 
more  good  stations  epicenter  location  MSE's  of  less  than 

200  km  are  obtainable.  Furthermore,,  if  a  region  can  be  cali¬ 
brated,  say,  by  a  number  of  shots  spaced  throughout  the  region, 
then  it  is  possible  to  obtain  location  MSE's  of  less  than  50  km2. 
In  this  case  there  is  essentially  no  source  bias  and  the  errors 
are  duo  to  either  measurement  errors  or  higher  order  model  in¬ 
adequacies  associated  with  the  stations. 


SUMMARY  OF  CONCLUSIONS 


The  following  paragraphs  provide  a  brief  summary  of  the 
most  important  conclusions  reached  in  this  report. 

1.  It  is  possible  to  develop  simplified,  practical, 
mathematical  models  for  the  prediction  of  teleseismic  travel 
times.  These  models  lead  to  a  non-trivial  definition  of  a 
"world-wide-average"  travel  time  function. 

2.  A  randomization  procedure  for  the  simplified  travel 
time  models  provides  the  basis  for  estimation  of  model  and 
location  parameters  under  two  extreme  conditions: 

(a)  All  errors  result  from  inaccuracies  of  measure¬ 
ment,  in  which  case  the  mean  square  time  error  is  minimized 
for  a  location  estimate,  and 

(b)  all  errors  result  from  model  inadequacies,  in 
which  case  the  mean  square  distance  error  is  minimized 
for  a  location  estimate.  In  practice,  errors  arise  from 
both  source,  and  either  estimation  technique  will  yield 
satisfactory  results.  Minimization  of  mean  square 
distance  errors  may  be  preferable  when  data  from  only 
very  good  stations  are  available. 

3.  When  estimating  travel  times  using  combined  data  from 
explosions  and  earthquakes,  all  the  sources  should  be  treated 
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in  the  same  manner.  That  is,  the  epicenters  should  be  esti¬ 
mated  for  the  explosions.  Utilization  of  known  location 
parameters  for  the  explosions  can,  in  this  case,  lead  to  bias 
in  the  travel  time  estimates. 

4.  The  1968  p-travel  times  (Herrin,  et  al. ,  1968)  are 
essentially  unbiased  and  statiscally  equivalent  to  the  most 
complete  set  of  shot  derived  travel  times  (Veith  and  Clawson, 
1972) . 

5.  The  estimated  station  corrections  given  by  Herrin 

and  Taggart  (1968a)  are  probably  adequate  for  those  cases  based 
on  50  or  more  observations. 

6.  Source  location  error  is  essentially  confounded  with 
the  first  order  azimuthal  source  terms.  This  effect  (source 
bias)  may  be  the  largest  cause  of  error  in  estimating  epicenters. 

7.  It  is  probably  useful  to  re-estimate  station  corrections 
to  remove  any  bias  due  to  mislocation  of  the  sources  and  to 
allow  for  distance  dependence. 

8.  With  "good"  station  corrections,  methods  can  probably 
be  developed  for  estimating  epicenters  with  mean-square-errors 
of  less  than  200  km^,  provided  the  network  consists  of  25  or 
more  well  distributed  stations.  These  methods  will  require  the 
use  of  deterministic  source  models  or  the  proportional  method 
for  approximating  the  source  interaction  term. 


TABLE  1* 


SUMMARY  OF  EXPLOSION  DATA 


Site  or  n  +■ 

No,  of 
Events 

No.  of 

Observations 

Total 

Teleseismic 

Semipalatinsk 

7 

528 

522 

Algeria 

6 

463 

442 

Novaya  Zemla 

9 

257 

251 

LONGS HOT 

1 

234 

232 

SHOAL 

1 

100 

46 

Nevada  Test  Site 

6 

531 

313 

SWORDFISH 

1 

48 

38 

CHASE  II 

1 

36 

33 

CHASE  III 

1 

53 

32 

CHASE  IV 

1 

48 

31 

CHASE  V 

i 

63 

31 

CHASE  VII 

1 

53 

33 

Bikini 

4 

21 

21 

Eniwetok 

_ 3 

19 

19 

Totals 

43 

2454 

2  044 

*From  Table  1,  Veith 

and  Clawson 

(1972). 

TABLE  2 


«iLALYSIS  OF  TUB  BXPTPSTnw  RESIDUALS  FOR  THE  HERR TT3  TRftVEI  TIMES 


Mean  values 

calculated 

per  5  dear, 

ee  blorlc 

Cell 

N 

Mean 

SD 

Variance 

SD  Mean 

T 

0-5 

100 

.262 

2.976 

8.856 

.298 

.880 

5-10 

135 

.026 

2.928 

8.573 

.252 

.103 

10-15 

104 

.955 

2.924 

8.551 

.287 

3.330 

15-20 

102 

1.05  0 

2.337 

5.462 

.231 

4.537 

20-25 

101 

.357 

1.722 

2.966 

.171 

2.086 

25-30 

100 

.083 

1.713 

2.933 

.171 

.485 

3  0-35 

109 

.367 

1.388 

1.927 

.133 

2.760 

35-40 

100 

.3  92 

1.237 

1.529 

.124 

3.170 

40-45 

86 

.112 

1.208 

1.459 

.13  0 

.857 

45-50 

72 

.072 

1.019 

1.039 

.120 

.601 

50-55 

56 

-.300 

.  988 

.976 

.13  2 

-2.272 

55-60 

70 

-.116 

1.085 

1.176 

.13  0 

-.893 

60-65 

84 

-.144 

1.212 

1.470 

.132 

-1.089 

65-70 

90 

-.  162 

1.091 

1.191 

.  115 

-1.410 

70-75 

80 

.041 

1.001 

1.003 

.112 

.368 

75-80 

77 

-.  296 

1.092 

1.192 

.124 

-2.380 

80-85 

85 

-.140 

.935 

.873 

.101 

-1.381 

85-90 

80 

-.043 

.  948 

.899 

.106 

-.401 

90-95 

68 

.068 

.  900 

.810 

.109 

.620 

95-100 

50 

.364 

.904 

.817 

.128 

2.847 

9l> 


TABLE  3 


PACIFIC -ATLANTIC  MODEL  STUDY  EVENT  REGION 


Event  Code 

Event  Reciion 

Amchitka 

Amchitka  Is.:  Aleutian  Trench 

Komandorskie 

Komandorskie  is . : 

Kuril 

Kuril  is.:  Kuril-Kamchatka  Trench 

Japan 

Japan:  Japan  Trench 

Ryukyu 

Ryukyu  Is.;  shoto  Trench 

Manila 

Philippines:  Philippine  Trench 

Western  New  Guinea 

Western  New  Guinea 

Tonga 

Tonga  Is.:  Kermadec  -  Tonga  Trench 

Chile 

Chile:  Peru-Chile  Trench 

Lima 

Peru:  Peru-Chile  Trench 

Columbia 

Columbia;  Pacific  Coast 

Guatemala 

Guatemala:  Pacific  coast 

Mazatlan 

Mexico:  Pacific  Coast 

Kod iak 

Kodiak  is. 

Unalaska 

Unalaska  Is.:  Aleutian  Trench 

Mariana 

Mariana  Is.:  Mariana  Trench 

Ogasawara 

Japan.*  Japan  Trench 

Java 

Java:  Java  Trench 

Sandwich  Is. 

So.  Sandwich  Is.:  South  Sandwich  Trench 

Puerto  Rico 

Puerto  Rico:  Puerto  Rico  Trench 

TABLE  4 _ 

PACIFIC -ATLANTIC  MODEL  STUDY  STATIONS 


Code 


Station 


Code 


Station 


AFI  Afiamalu,  Samoa 

ATH  Athens,  Greece 

BAG  Baguio,  Philippines 

BRK  Berkley,  California 

BOG  Bogota  (WWNSS),  Columbia 

BOK  Bokar,  India 

BUL  Bulawayo,  Southern  Rhodesia 
CAN  Canberra,  Australia 
CTA  Charters  Tower,  Australia 
CMO  College,  Alaska 
CSC  Columbia,  So.  Carolina 
DAL  Dallas,  Texas 

GUA  Guam,  Guam 

HLW  Helwan,  Egypt 

HVO  Hawaiian  Volcano  Obs., 
Hawaii 

KIS  Kishinev,  USSR 

LAN  Lanchow,  China 

LPS  LaPalma ,  El  Salvador 

LAR  Laramie,  Wyoming 

LHA  Lhasa,  China 

LPB  La  Paz,  Bolivia  (WWSS) 

LWI  Lwiro ,  Congo 

MED  Medan,  Sumatra 

MNT  Montreal,  Canada 
MOS  Moscow,  USSR 

MBC  Mould  Bay,  Canada 


PAL  Palisades,  New  York 

PER  Perth,  Australia 

PVC  Port  Vila,  New  Hebrides 

QUE  Quetta,  W.  Pakistan 

REY  Reykjavik,  Iceland 

ROM  Rome,  Italy 

SAP  Sapporo,  japan 

SEM  Semipalatinsk,  USSR 

SEN  Sendai,  Japan 

SHI  Shiraz,  Iran 

SIM  Simferopol,  USSR 

SOD  Sodankyla,  Finland 

SAN  Santiago,  Chile 

THU  Thule,  Greenland 

TIK  Tiksi,  USSR 

TOL  Toledo,  Spain 

TRN  Trinidid,  West  Indies 

UPP  Uppsala,  Sweden 

VAN  Vannovskaya,  USSR 

VCM  Vera  Cruz,  Mexico 

VIE  Vienna,  Austria 

VLA  Vladivostok,  USSR 

WAR  Warsaw,  Poland 

WEL  Wellington,  New  Zealand 

WIL  Wilkes,  Antarctica 

WIN  Windhoek,  South  Africa 


TABLE  5 


PACIFIC-ATLANTIC  MODET,  STUDY 


Cell 

Mean 

SD  Mean 

N 

20-25 

.  035 

.449 

26 

25-30 

.173 

.314 

20 

30-35 

.185 

.227 

25 

35-40 

.041 

.  193 

29 

40-45 

.129 

.162 

23 

45-50 

.018 

.147 

32 

50-55 

.016 

.106 

38 

55-60 

.040 

.078 

39 

60-65 

.008 

.063 

40 

65-70 

.006 

.066 

49 

70-75 

.006 

.119 

53 

75-80 

-.021 

.172 

51 

80-85 

-.049 

.177 

56 

85-90 

-.121 

.205 

57 

90-95 

-.073 

.256 

44 

95-100 

-.  060 

.  256 

51 

100-105 

.004 

.259 

51 

TABLE  6 


STATIONS  EMPLOYED  IN  THE  AMCHITKA  MONTE  CARLO  STUDY 


Code  Station 


Code  Station 


API 

ATH 

BAG 

BRK 

BOG 

BOK 

CAN 

CTA 

CMO 

CSC 

DAL 

GUA 

HLW 

HVO 

KIS 

IAN 

LPS 

LAR 

LHA 

MED 

MNT 

MOS 

MBC 


Afiamalu,  Samoa 
Athens,  Greece 
Baguio,  Philippines 
Berkely,  California 
Bogota  (WWNSS ) ,  Columbia 
Bo* ,  India 
Canberra,  Australia 
Charters  Tower,  Australia 
College,  Alaska 
Columbia,  South  Carolina 
Dallas,  Texas 
Guam,  Guam 
Helwan,  Egypt 

Hawaiian  Volcano  Obs., Hawaii 

Kishinev,  USSR 

Lanchow,  China 

LaPalma,  El  Salvador 

Laramie ,  Wyoming 

Lhasa ,  China 

Medan,  Sumatra 

Montreal,  Canada 

Moscow,  USSR 

Mould  Bay,  Canada 


PAL 

PER 

PVC 

QUE 

REY 

ROM 

SAP 

SEM 

SEN 

SHI 

SIM 

SOD 

THU 

TIK 

TOL 

TRN 

UPP 

VAN 

VCM 

VIE 

VLA 

WAR 

WEL 


Palisades,  New  York 
Perth,  Australia 
Port  Vila,  New  Hebrides 
Quetta,  'West  Pakistan 
Reykjavik,  Iceland 
Rome,  Italy 
Sapporo,  japan 
Semipalatinsk,  USSR 
Sendai,  Japan 
Shiraz,  Iran 
Simferopol,  USSR 
Sodankyla,  Finland 
Thule,  Greenland 
Tiksi,  USSR 
Toledo,  Spain 
Trinidad,  west  Indies 
Uppsala,  Sweden 
Vannovskaya,  USSR 
Vera  Cruz,  Mexico 
Vienna,  Austria 
Vladivostok,  USSR 
Warsaw,  Poland 
Wellington,  New  Zealand 


TABLE  7 


AMCHITKA  MONTE  CARLO  STUDY  RESULTS 


n 

6 

N-S  Mean 

N-S 

&  N-S 

E-W  Mean 

^  E-W 

?E-w 

Source  Estimated* 

128 

1.0 

-4.261 

17. 93 

17.93 

-0.560 

17.16 

17.16 

48 

0.8 

-3.876 

12.11 

15.14 

1.810 

16.28 

20.35 

48 

0.08 

0.126 

1.269 

15.86 

-0.015 

1.763 

16.60 

Source  Not  Estimated* 

88 

1.0 

-0.715 

3.988 

3.988 

-0.179 

3.431 

3.43. 

♦North  and  east  are  plus  while  south  and  west  are  minus. 


TABLE  8 


TEST  NETWORKS  FOR  AMCMITKA  MONTE  CARLO  STUDY 


Network 


Station 


I  Berkeley,  California 

Hawaiian  Volcano  Obs.,  Hawaii 

Laramie,  Wyoming 

Mould  Bay,  Canada 

Sapporo,  Japan 

Sendai,  Japan 

Thule,  Greenland 

Tiksi,  USSR 

Vladivastok,  USSR 


II  Afiamaul,  Samoa 

Columbia,  South  Carolina 
Lhasa,  China 
Moscow,  USSR 
Palisades,  New  York 
Port  Vila,  New  Hebrides 
Uppsala,  Sweden 
Vera  Cruz,  Mexico 


III  Berkeley,  California 

Laramie ,  Wymoming 
Mould  Bay,  Canada 
Thule,  Greenland 

_ Tiksi.* -USSR _ 

IV  Hawaiian  Volcano  Obs.,  Hawaii 

Sappor  o,  Japan 

Sendai,  japan 
Tiksi,  USSR 
Vladivastok,  USSR 


Oj(p 


TABLE  9 


AMCHITKA  TEST  NETWORK  ERRORS 


Network 

No.  of 
Stations 

Source  not 
Estimated* 

Source 

Estimated* 

I 

9 

17.21 

-0.95 

-0.15 

0.05 

II 

8 

23.15 

-1.34 

-1.31 

0.74 

III 

5 

18.05 

-0.06 

0.45 

-0.06 

IV 

5 

16.61 

-1.43 

72.93 

-3.11 

♦North  and 

east  are  plus 

while  south  and 

west  are  minus, 

c\1 


Network 

V 


VI 


TABLE  10 

LONGS MOT  STUDY  NETWORKS 


Code 

Station 

Code 

ALQ 

Albuquerque,  New  Mex. 

LAR 

BNS 

Bensberg,  Germany 

MAT 

CAN 

Canberra,  Australia 

MNT 

CTA 

Charters  Tower,  Australia 

NUR 

CMC 

Copper  Mine,  Canada 

PAS 

COL 

College  outpost,  Alaska 

PMG 

CPO 

Cumberland  plateau  S.O., 

Tennessee 

PRU 

DAL 

Dallas,  Texas 

QUE 

DH" 

Delhi,  New  York 

RES 

FBC 

Frobisher,  Canada 

SHI 

ISO 

Isola,  France 

SOD 

JER 

Jerusalem,  Israel 

STR 

KJN 

Kajaani,  Finland 

TOO 

KHC 

Kasperske  Hory, 

TRO 

Czechoslovakia 

UPP 

KEV 

Kevo,  Finlan 

VAL 

KIR 

Kiruna,  Sweden 

WMO 

KON 

Kongsberg,  Norway 

BKS 

Byerly,  California 

LON 

BUT 

Butte,  Montana 

MAT 

COL 

College  outpost,  Alaska 

MN~ 

EUR 

Eureka,  Nevada 

PRS 

FBC 

Frobisher,  Canada 

PAS 

GOL 

Golden,  Colorado 

SOD 

KEV 

Kevo,  Finland 

SPO 

KIR 

Kiruna,  Sweden 

TFO 

LAR 

Laramie,  Wyoming 

LAW 

Lawrence,  Kansas 

TUC 

LC" 

Las  Cruces,  New  Mex. 

UBO 

MHC 

Lick  (Mt.  Hamilton), 

California 

VIC 

WMO 

Station 

La  ramie ,  Wyoming 
Matsushiro,  Japan 
Montreal,  Canada 
Nurmi jarvi,  Finland 
Pasadena,  Calif, 
pt.  Moresby,  New 
Guinea 
pruhonice,  Czech. 
Quetta,  W.  Pakistan 
Resolute  Bay,  Cana. 
Shiraz,  Iran 
Sodankyla,  Finland 
Strasbourg,  France 
Toolangi,  Australia 
Tromsoe ,  Norway 
Uppsala,  Sweden 
Valentia,  Ireland 
Wichita  Mts.  Seis. 
Obs.,  Oklahoma 


Longmire,  Wash. 
Matsushiro,  Japan 
Mina,  Vevada 
Paraiso,  Calif. 
Pasadena,  Calif. 
Sodankyla,  Finland 
Spokane,  Washington 
Tonto  Forest  Seis. 
Obs.,  Arizona 
Tucson,  Arizona 
Uinta  Basin  Seis. 
Obs . ,  Utah 
Victoria,  Canada 
Wichita  Mts.  Seis. 
Obs.,  Oklahoma 


aablk 


LONGS HOT  STUDY  NETWORK  ERRORS 


Network 

North-South 

Error* 

East-West 

Error* 

Source 

Not  Es¬ 
timated 

No  Station 
Corrections 

19.86 

-1.66 

V 

Station 

Corrections 

21.08 

-1.5C 

Source 

Estimated 

No  Station 
Corrections 

“19.74 

1.28 

Station 

Corrections 

-35.51 

10.19 

Source 

Not  Es¬ 
timated 

No  Station 
Corrections 

30.22 

-2.57 

VI 

Station 

Corrections 

21.03 

2.05 

Source 

Estimated 

No  Station 
Corrections 

-9.40 

-53.11 

Station 

Corrections 

-66.28 

6.72 

♦North  and  east  are  plus  while  south  and  west  are  minus 


(spuosw)  ivnaisau 

too 
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Figure  3  AMCHITKA  MODEL 
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Figure  5  COMPUTED  RESIDUALS 


Figure  6  MODEL  RESIDUALS  VS.  AZIMUTH 
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